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THE THERMODYNAMIC BASIS FOR THE 
MINERAL FACIES CONCEPT 


JAMES B. THOMPSON, Jr. 


ABSTRACT. For a given petrologic system, the curves of univariant equilibrium on a 
P, T diagram may be regarded as a “petrogenetic grid” (Bowen, 1940) defining the 
mineral facies for that system. 


If only crystalline solids are involved it can be shown that most curves of this sort 
are essentially straight lines through any geologically probable range of temperature and 
pressure. 


Mineralogic reactions involving the so-called “volatile” components are herein treated 
by considering the system as “open” to these components. Four limiting cases are dis- 
cussed and their geologic applications reviewed briefly. The limiting cases are: 


1. The “volatile” component in question is present as a pure phase at the same 
pressure and temperature as the surrounding rock. This assumption yields a maximum 
equilibrium temperature for a given pressure, and is roughly equivalent to considering 
the system as closed with the volatile component present in excess. 


2. The rock is in equilibrium with a fissure system containing a dilute aqueous 
solution. This method is most applicable when the “volatile” in question is water itself, 
and indicates, in general, a decrease in dehydration temperature with depth, 


3. The gradient of chemical potential of the “volatile” components has reached a 
value such that there is no longer any tendency for their vertical diffusion. 


4. The chemical potential of the “volatile” component is sufficiently low that no 
solid phase containing it is stable. This is analogous to a closed system devoid of the 
volatile component. 


NOTATION 
any extensive property 
internal energy 
enthalpy 
Gibbs free energy 
entropy 
volume 
1 absolute temperature 
P pressure 
C, heat capacity at constant pressure 
a thermal expansion 
B compressibility (isothermal) 
p density 
Subscripts 
S solid phases 
F fluid phases 
a, b, c, individual phases 
i, j, k, individual components 
Nia number of moles or gram-formula units of component i in phase a 
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Xia partial molar X of component i in phase a 

Nia mole fraction of component i in phase a 

AX X of products less X of reactants 

Mia (—G,,) chemical potential of component i in phase a 
g acceleration of gravity 

R gas constant 


INTRODUCTION 


There is considerable textural and mineralogic evidence (Eskola, 1915, 
1920; Tilley, 1924; Turner, 1948) that many rocks represent a state of 
chemical equilibrium, or at least a close approximation thereto, as of their 
time of formation. Where this is true, the mineralogy of a given rock should 
be the key, within limits, to the conditions prevailing when it formed. This 
is the basis of Eskola’s concept (1920) of mineral facies. Bowen (1940) has 
suggested that a pressure-temperature diagram upon which have been plotted 
the curves of univariant equilibrium appropriate to a given bulk composition 
rock might serve as a “petrogenetic grid.” The areas between the lines of 
univariant equilibrium would then be the fields of stability of the various 
mineral facies for that bulk composition. With this end in view experimental 
investigation of the field of stability of \ arious mineral species has been under 


way for several years (Bowen and Tuttle, 1949; Yoder, 1950a, 1950b, 1951, 
1952) : 


A simple P. T. grid of this sort, however, presupposes that the bulk com- 
position of the rock may be regarded as fixed, an assumption that cannot be 
made safely with regard to highly mobile components such as water and 
carbon dioxide. Petrographic evidence indicates that these components are 
readily gained or lost in petrogenetic processes. Thermodynamically, this 
means that we must regard a rock or petrologic “system” as open to certain 
mobile components, notably water and carbon dioxide. The field of stability 
of a mineral containing a mobile component is then dependent not only upon 
temperature and pressure but also upon the chemical potential (partial molar 
free energy) of the mobile component in the immediate environment. In a 
completely closed system the chemical potential of any component is de- 
termined by the pressure, temperature, and bulk composition, If a system is 
open to a certain component, however, its chemical potential is independent of 
temperature and pressure, owing to the possibility of transfer of the mobile 
component between a given system and its environment, The field of stability 
of a mineral facies in a rock containing m mobile components is thus, in 
general, (m + 2)-variant, that is, dependent on pressure, temperature and 
the chemical potential of the mobile components. The boundary between two 
mineral facies is then no longer a univariant curve on a P-T diagram but 
an (m + 1)-variant surface or hypersurface in (m + 2)-dimensional space. 

Fortunately, few metamorphic changes involve gain or loss of more than 
one mobile component at a time. Commonly the only mobile component 
involved is water or carbon dioxide. Such equilibria are divariant surfaces 
in a 3-dimensional diagram, where pressure, temperature, and the chemical 
potential of the mobile component are the coordinates. If, in any instance, the 
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chemical potential gradient of the mobile component may be regarded as 
dependent on the temperature and pressure distribution in the Earth’s crust, 
the equilibrium may be regarded, to the extent that this is true, as univariant 
and dependent only on pressure and temperature. A possible mechanism by 
which this last might come about is through attainment of osmotic equilibrium 
with a hydrothermal fissure system. Another possibility is that free vertical 
diffusion of the mobile components has taken place so that an equilibrium 
chemical potential gradient has been attained for each. Both possibilities will 
be considered in the following pages, as well as certain features of univariant 
equilibria in closed systems. 

Earlier work.—The work of Johnston and Niggli (1913) and V. M. 
Goldschmidt (1912) is notable among early investigations of the physico- 
chemical principles involved in metamorphic processes, A revival of interest 
in such matters has been indicated by the chapters on chemical thermody- 
namics in recent texts on petrology and geochemistry, chiefly those by Turner 
and Verhoogen (1951), Mason (1952), Barth (1952), and Ramberg (1952). 
Certain important aspects of phase equilibrium in open systems have been 
discussed by Korzhinsky (1936, 1950). The calcite-wollastonite equilibrium 
has been restudied by Danielsson (1951) from the point of view of osmotic 
equilibrium, and Ramberg (1944, 1952) has considered free vertical diffusion 
of mobile components as a factor in metamorphic processes. 

Recent experiments in hydrothermal systems, notably by Bowen and 
Tuttle (1949), Roy and Osborne (1950), and Yoder (1950a, 1950b, 1952), 
has also provided much-needed data and stimulated interest in the réle of 
water in petrologic processes. Yoder and Weir (1951) have considered the 
effect of pressure on metamorphic processes, and Yoder (1952) has con- 
sidered the fields of stability of certain hydrous minerals from the point of 
view of closed systems with “excess” or “deficient” water. 

Readers familiar with the above literature are probably aware of funda- 
mental differences in the approaches of the various authors, and of the difh- 
culties involved in correlating field observations with either experimental 
data or theoretical discussions. It is the hope of the present author to be able, 
in this paper, to resolve some of the inconsistencies in the current literature 
and to narrow somewhat the existing gaps between isograds, P-T curves, and 
differential equations. 


UNIVARIANT EQUILIBRIUM 
Solid phases of fixed composition—The equation of univariant equili- 
brium for reactions at constant pressure and temperature may be written 
AG = 0 = AE + PAV - TAS (1) 
or AG = 0 = AH- TAS (la) 
The reaction is written so that AS is positive. From the first and second laws 
of thermodynamics we have, for a reversible process 
d(AE) = Td(AS) — Pd(AV) 


hence, maintaining equilibrium, we have 


3 
dT AG—0 AV 
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or, alternatively, 


( —) AH 
aT ~ TaVv 
AG=0 : 
a relationship commonly known as the Clapeyron equation. 
For most solid-solid reactions at moderate (room temperature) and high 
temperatures it can be shown that equation (1) is, for most practical purposes, 
the equation of a straight line on a P-T diagram, the slope of this line, for a 
given reaction, being given by (3). 
The volume of most solids varies little with temperature and pressure. 
From available data (Birch et al., 1942, p. 30-37) we find that, at 25°C. and 
one atmosphere, the volume expansion, a, defined by the relation 


en. (5) 
V \ oT P (4), 


is between 1 x 10° per degree and 4 x 10° per degree for most rock-forming 
minerals, and that the isothermal compressibility, 8, defined by the relation 


= 
a (ae T (5), 


is between 0.5 x 10°° per bar and 3 x 10° per bar. The variation of « and 
8 with temperature and pressure has been discussed at some length by Birch 
(1952), both quantities, in general, increasing slightly with rising tempera- 
ture and decreasing slightly with rising pressure. Nevertheless, these varia- 
tions, as well as the initial values of « and 8, are so small that we may neglect 
them for most “chemical” purposes despite their importance in geophysical 
calculations. Since most solids have about the same thermal expansion and 
compressibility, it is evident that we are even safer in assuming constancy 
of AV for any univariant equilibrium than of V for any of the individual 
participating phases. The volume changes of the reactants, in short, are largely 
cancelled by corresponding volume changes on the part of the products, so 
that AV is virtually constant. This constancy has, in fact, been effectively 
demonstrated by Yoder and Weir (1951) for the specific reaction: 2 jadeite 

- nepheline + albite.t Available compressibility and thermal expansion data 
indicate that it is probably a safe assumption for any univariant equilibrium 
involving only solid phases. 

Molar volumes may be determined readily and with considerable ac- 
curacy, when necessary, from either specific gravity or unit cell measurements. 
The minerals in table 1 are listed in order of increasing volume per gm-atom. 
The volumes, per gram-atom, of most mineral oxides, oxy-salts, and silicates 
lie in the range 4.5 to 9.0 ce (table 1). Values of 7.5 cc/gm-atom or more 
are found chiefly in silicates with tetrahedral frameworks (feldspars, feldspa- 
thoids, silica minerals, and cordierite), and values of 5.0 cc/gm-atom or less 


* Although demonstrating the virtual constancy of AV, Yoder and Weir draw therefrom 
an erroneous conclusion. Specifically, the constancy of AV, and hence of the effect of 
pressure on AG, does not mean that the effect itself is zero or negligible. AV, for the 
reaction: 2 jadeite = nepheline + albite, although virtually constant, is neither zero 
nor negligible. 


(3a) 





Basis for the Mineral Facies Concept 69 


occur chiefly in hydroxides. Volume changes accompanying solid-solid transi- 
lions are commonly on the order of 0.5 to 1.0 cc/gm-atom, Extreme examples 
are the reaction diamond — graphite, and the reaction 2 jadeite — albite + 
nepheline. Among polymorphous pairs, that form which lies highest on the 
list in table 1, is the one favored by an increase in pressure. The greater the 
volume-difference, the more marked the effect of a given increase in pressure. 

























































































TABLE 1 
Volumes (25°C.) 
Molar Gram-atomic 
Mineral Formula Volume (cc) Volume (cc) 
Diamond c 3.4 3.4 
Diaspore HAIO. 17.6 4.4 
Gibbsite A1(OH)s 32 4.6 
Brucite Mg(OH)s 24.2 4.8 
Corundum Al,Os 26 LY 
Graphite Cc 5.3 5.3 
Staurolite Fe(OH) o* Ali (Sig) 20a ...esceeeeeeeeeet 111 5.3 
Kyanite Al.(Si0,)0 44 5.5 
Lawsonite CaAl.(Sis0:) (OH)2*H2O «0.0.0.0... 104 5.5 
Chloritoid Te(OH) 2*Ale(SiOg)O  .....ceceserereeee 72.5 5.6 
Periclase MgO 11.2 5.6 
Pyrope MgsAle(SiO,.)s 113 5.7 
Almandine FesAle(SiO,)s 115 5.8 
Kaolinite Al. (SizOs) (OH), eecccccccccceocoscescoscooss 99 5.8 
Chlorite MgsAl( AISisOw) (OH) 5 .....-0eceese+2 211 5.9 
Jadeite NaAl(SisOs) 60 6.0 
Serpentine Mige(SisOs) (OH) 4 -n-cccccscecsscccsescocees 109 6.1 
Sillimanite Al(AISiO;) 44 6.1 
Forsterite Mge(SiO,) 44 6.3 
Enstatite Mg(SiOs) 31.5 6.3 
Grossularite CasAle(SiO,)s 126 6.3 
Clinozoisite Casble(SiQs) (OH)  ..cccrcceccccossesscsees 141 6.4 
Pyrophyllite Ale(SigO10) (OH). eneuseeensvencetecenecsnss 128 6.4 
Tale Mega(SisOro) (OH) 2 cecsscscecssoscersvesees 136 6.5 
Andalusite AlAI(Si0,)O 52 6.5 
Muscovite KAle( AISisO1w) (OH) 2 ..ecccceseeceeeeee 142 68 
Phlogopite KMgs( AlSisO10) (OH) 2 ...secesceeeeeeee 150 6.8 
Aragonite Ca(COs) 34 6.8 
Calcite Ca(COs) 37 7.4 
Quartz (Low) SiO, 22.4 7.5 
Analcite | ©. Sen enee 98 7.5 
Gehlenite Caz( Al.SiO;) 91 7.6 
Anorthite Ca(AbLSi:0s) 101 7.8 
Albite Na(AISisOs) . 102 7.8 
Wollastonite Ca(SiOs) 39.5 7.9 
Nepheline Na(AISi0,) 56.5 8.1 
Orthoclase K( AlSisOs) 105 8.1 
Cordierite Mge( ALSis0;s) 236 8.1 
Tridymite (Low) SiO: 25.5 8.5 
Cristobalite (Low) SiO. 26.0 8.7 











Volumes for garnets are based on unit-cell dimensions given by Fleischer (1937), all 
others are based on unit-cell volumes listed by Fairbairn (1943). 
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The molar and gram-atomic entropies (at 25°C. and 1 atm.) of some 
minerals have been listed in table 2. The listing is in order of increasing en- 
tropy per gram-atom. The relation of entropy to temperature is analogous 
(but opposite in sign) to that of volume to pressure, Among polymorphous 
pairs, that form which is lowest on the list is the one most favored by an 
increase in temperature. The magnitude of this effect is proportional to the 
entropy-difference. At room temperature and atmospheric pressure the gram- 
atomic entropies of most rock-forming minerals are between 100 and 225 
decijoules. 

Materials of low gram-atomic entropy are commonly materials of low 
gram-atomic volume, as is evident from comparison of tables 1 and 2. There 
are exceptions, an important one being andalusite-sillimanite, but these ex- 
ceptions are limited chiefly to cases where the volume and entropy changes 
are relatively slight. This means that an increase in pressure generally has 
an effect, on a solid-solid transition, opposite to that produced by an increase 
in temperature. The point should be emphasized because it is contrary to the 
impression given in many discussions of metamorphic processes. Most equili- 
brium curves for solid-solid reactions will thus have positive slopes on a P-T 
diagram, The exceptions, in fact, are sufficiently rare that it is fairly safe 
to assume, until shown otherwise, that the phase or “facies” with the smaller 
volume will also have a smaller entropy and hence be the “low-temperature” 
form. This is particularly so if the volume change per gram-atom is relatively 
large, and is true in all instances known to the writer where AV exceeds 0.5 
cc/gm-atom. 

TABLE 2 
Entropies (25°C.) 

Molar Entropy 
(dj/deg) 


24.4+0.2 
56.9+0.2 
701+4.0 
523 +6.0 
838 + 3.0 
932+4.0 
961+4.0 
631+2.0 
1331+8.0 
679+4.0 
952 +8.0 
274+6.0 
418+4.0 
426+ 4.0 
433+8.0 


Gram-atomic 
Entropy (dj/deg) 


24.4+0.2 
56.9 +0.2 
100.1 +0.6 
104.6 + 1.2 
104.7+0.4 
116.5+0.5 
120.1+0.5 
126.2+0.4 
133.1+0.8 
135.8+0.8 
136.0+ 1.1 
137+3.0 
139+1.3 
142+1.3 
144+2.7 


Mineral Formula 








Diamond 
Graphite Cc 
Gibbsite Al(OH): 
Corundum 

Kyanite 

Andalusite 
Sillimanite 

Brucite 

Jadeite 

Enstatite 

Forsterite 

Periclase 

Quartz (Low) 
Cristobalite (Low) * 
Tridymite (Low) * 
Portlandite 


(AlSiOs) 
Mg(OH). 
NaAl (SizOs) 


Albite 
Wollastonite 
Nepheline 
Aragonite 
Calcite 

Lime 

Water (liq.) 


( ‘a(OH )e 
Na( Al1Si,05 ) 
CaSiO, 
Na(AISiO,) 
Ca(COs) 





730 + 40.0 
2059 + 13.0 
820+ 8.0 
1218+8.0 
887 + 13.0 
929+8.0 
410+8.0 
701 + 13.0 


146+8.0 
158.4+1.0 
164+ 1.6 
174+1.1 
177+2.6 
186 + 1.6 
205 + 4.0 
234+ 4.0 





Entropy of Nepheline from Kracek et al. (1951) 
Data converted to decijoules/deg (1 dj = 41.84 cal.). 
* Relative positions reversed for “high” forms. 


, all other data from Kelley (1950). 
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AS for most solid-solid transitions is commonly between 5 and 10 deci- 
joules per gram-atom per degree. Extreme examples, again, are 32.5 deci- 
joules/gm-atom/degree for diamond — graphite and 30.8 decijoules/gm- 
atom/degree for 2 jadeite — nepheline + albite. The thermochemical data 
(table 2) indicates that the uncertainty in AS may be from + 1 to + 5 
decijoules/gm-atom/degree. This should be borne in mind in all instances 
where AS /gm-atom is relatively small. 

The entropy of a crystalline solid, as given in thermochemical tables, is 
based upon low-temperature heat capacity measurements and, strictly speak- 
ing, is an entropy difference; it is the entropy of the crystal at 25°C, less its 
entropy at absolute zero. It is commonly assumed, for chemical purposes, 
that the latter is zero, but this is not necessarily so. Although the vibrational 
entropy of a crystal drops to zero as the temperature approaches absolute 
zero, other contributions to the entropy (as, for example, the disordering of 
aluminum and silicon, in high-temperature feldspars) may readily be “frozen 
in.” If this has happened, and is not taken into account, the heat capacity 
measurements will give too low a value by that part of the total entropy which 
had been “frozen in.” Similarly, heat capacity measurements on a mix-crystal 
such as an intermediate olivine will neglect the entropy of mixing. Fortunately 
it is pessible in most instances to anticipate difficulties of this sort from a 
knowledge of the crystal structure, and to make a reasonable estimate of the 
entropy-difference involved, though this is beyond the scope of the present 
paper. 

The effect of pressure on entropy at constant temperature is given by 


the relation 
——- == os —_ —_— — a 
oP T oT P (6) 


Since the thermal expansion of solids is negligible for our purposes, it is clear 

that pressure has little effect on S, and even less on the AS of a reaction in- 

asmuch as the thermal expansions of reactants and products will be similar. 
The variation of entropy with temperature is given by the relation 


( a 1 (5) Cp 

oT P T \ oT p T (7) 
where C, is the heat capacity at constant pressure. At room temperature and 
one atmosphere the heat capacity of most rock-forming minerals is between 
150 and 210 decijoules/gm-atom/degree. At lower temperatures heat capaci- 
ties decrease rapidly, approaching zero as the temperature approaches absolute 
zero. At higher temperatures the heat capacities of virtually all solids ap- 
proach, asymptotically, values of about 250 decijoules/gm-atom/degree, In 
other words, the gram-atomic heat capacities of solids approach, at high 
temperatures, the value of 3R (249.42 decijoules/gm-atom/degree) indicated 
by the empirical law of Dulong and Petit (Slater, 1939, p. 213-214) as well 


as by the later theories of Einstein and Debye (Slater, p. 222-255). Strictly 
speaking, the Einstein and Debye theories predict that the heat capacity at 
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constant volume, C,, rather than C,, should approach 3R/gm-atom, but as 
we can see from the relation 


a 


. a 
Cp — Cy + Ive (8) 


these do not differ greatly in solids. The gram-atomic heat capacities of most 
rock-forming minerals approach 3R in the temperature range from 600 to 


1000°C. Since 
(=) _ ACp 
oT} p ‘ (9) 


it is clear that AS for an isochemical solid-solid transition may be regarded 
as constant at temperatures above 1000°C. Actually this is sufficiently true 
even at room temperature inasmuch as the deviations from 3R are about the 
same for any two polymorphs or isochemical aggregates, and T is a relatively 
ACp 
T 
joules /gm-atom/degree* and becomes rapidly less with rising T, owing to both 
the increase in T itself and the fact that ACp approaches zero. We could take 
this into account in the equations of univariant equilibrium, but in view of 
the initial uncertainty in determining AS there is little point in so doing. 

It is evident, then, that if we restrict the discussion to temperatures above 
25°C. (safe enough geologically) we may assume, as a first approximation, 
that both AS and AV are constant for solid-solid reactions, despite the fact that 
the entropy of any one of the participating phases is dependent on tempera- 
ture. 


large number. For most reactions that concern us is less than .05 deci- 


It can also be shown that the ratio AS/AV for a solid-solid reaction mav 
also be regarded as constant, since, 





AV >. ae 


ee AS 
r) (=) | _ 4S [ “a A(Va) 


oP 


_ AAV] £4 34 ae 
oT 


a _/P (11) 


For most reactions, the bracketed terms on the right hand side of each 
of these relations tend to cancel each other because entropy and volume vary 
sympathetically. The largest of these bracketed terms, at least at low tempera- 
tures, will generally be the heat-capacity term in equation (11). This means 
that we may anticipate that some solid-solid equilibrium curves will show, at 
relatively low temperatures, a slight concavity toward the pressure axis, since 
the high-entropy phase has, in general, the higher heat capacity. 

From equation (1) and equation (3) it is evident that constancy of AS 


@(=) _ aS[7ACp A(Va) 
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and AV implies constancy of AE, at least along the equilibrium curve. This 
is also evident from the relations 


dE Saale 
5P — PVB —~ TVe 


ee 
a Cg 
oT] p (13) 


From equations (7) and (9) it is further evident that AH may be regarded 
as independent of temperature at constant pressure. This suggests that if AH 
and AS at 25°C, and one atmosphere are known we might obtain the transition 
temperature at one atmosphere by means of equation (la). Unfortunately, 
however, relative enthalpies, as determined from heats of solution, are gen- 
erally not of sufficient accuracy to make such a calculation worthwhile, and 
equation (la) is useful chiefly as a way to determine the latent heat of a 
reaction when the entropy change and transition temperature are known, 
at least until more accurately determined heats of solution are forthcoming. 
Most measurements of heats of solution of minerals are given as good to about 
the nearest hundred calories per gram-atom, that is, to within +50 calories 
(2000 decijoules) per gram-atom. AH is then good to within +4000 deci- 
joules per gram-atom. If AS is 10 decijoules/gm-atom/degree (30 decijoules/ 
gm-atom/degree is an exceptionally large value) this means an uncertainty 
of +400°C. in determining the transition temperature by equation (la), 
neglecting any error in AS. For an exceptionally accurate heat of solution, 
or an exceptionally high AS, or both, this uncertainty is reduced. With present 
data, then, it is often better to estimate the position of an equilibrium curve 
by an “educated geological guess” than to employ equation (la), although, 
since the errors in AS/AV are considerably less, we may estimate its slope 
with some accuracy by means of equation (3). 

Non-first-order transitions.—Certain transitions in solids, notably those 
of the order-disorder type, may take place with no discontinuity in the first- 
order derivatives of the free energy, that is, without discontinuities in entropy 
and volume, though these properties may change rapidly over a narrow range 
of temperature or pressure. In such an instance AS and AV are both zero, 
and the Clapeyron equation (3) is indeterminate. When discontinuities ap- 
pear in second- or higher-order derivatives of the free energy, the pressure 
dependence of the transition temperature may be determined from equations 
derived by Ehrenfest (Guggenheim, 1950, p. 775), but practical difficulties 
arise in measuring such discontinuities. 

An alternate approach is to consider the rapid changes in entropy and 
volume as taking place abruptly rather than over an interval. This is, in 
effect, considering a hypothetical equilibrium between (in the case of an order- 
disorder transition) the disordered state and a partially ordered state. Both 
phases in such an equilibrium are metastable or even unstable relative to 
some state with an intermediate degree of order, but this is of no consequence 
for our purposes. Such an approximation is useful because the pressure-de- 
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pendence of such a metastable equilibrium can be evaluated, and because the 
P-T curve for the metastable or unstable transition must always be within the 
band representing the actual transition interval. The value of this method is 
greatest for relatively abrupt transitions, specifically those for which the 
transition band on a P-T diagram is relatively narrow. 

Phases of variable composition—An analagous situation occurs when, 
with phases of variable composition, a given mineralogic change takes place 
over an interval rather than at a fixed P or T. Examples are common and are 
to be anticipated wherever an intermediate member of a solid solution series 
or a liquid phase is involved, and where the total number of phases involved 
is less than the maximum permissible for univariant equilibrium. Thus the 
reaction: 

orthopyroxene — clinopyroxene 

takes place over a temperature interval, for intermediate Fe-Mg ratios (Bowen 
and Schairer, 1935), but the reaction: 

olivine + orthopyroxene + tridymite — olivine + 

clinopyroxene + tridymite 
is abrupt at constant pressure owing to the phase-rule restrictions, In the 
first reaction the transition interval for any given composition would be a 
band on a P-T diagram analogous to that for order-disorder transitions. 

A mineral aggregate, on passing through such a transition interval, with 
maintenance of equilibrium, undergoes two second-order transitions, one on 
entering the interval and one on leaving it, but no first-order transition, The 
discontinuities in heat capacity and thermal expansion may be used as in 
other second-order transitions, but it is more convenient, for most purposes, 
to consider a related metastable or unstable equilibrium. In the above example 
it would be the P-T curve along which an orthopyroxene has the same molar 
free energy as a clinopyroxene of identical composition. This can be no more 
than a partial equilibrium since the chemical potentials of the pyroxene end- 
members must be different in the two phases, but this does not prevent its being 
defined thermodynamically. The partial equilibrium always lies within the 
actual transition interval. 

In this paper the equations for most reactions will be written in terms 
of pure minerai end members. An actual reaction involving mix crystals will 
in general behave differently owing to the effects of mixing upon the thermo- 
dynamic properties. In some minerals, notably in the Fe-Mg olivines and 
orthopyroxenes (Sahama and Torgeson, 1949) and in the plagioclase feld- 
spars (Bowen, 1928, p. 179), there is evidence that the solid solutions are 
ideal or nearly ideal. In an ideal mixture, a, the free energy of the mixture 
may be given in terms of the chemical potentials (y;°) of the pure end mem- 
bers by the equation 

G= xi Niatia — Xi Nia (pi” + RTinN j,) 
(14) 
It is probably not safe to assume, however, that many solid solutions are ideal 
mixtures. 
The effect of non-hydrostatic stress—For most chemical purposes equa- 
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tions (1) and (2) are adequate starting points. Both carry the assumption, 
however, that the work done in a chemical process is solely that involved 
in changes of volume against hydrostatic pressures. This cannot always be 
justified in a petrologic system inasmuch as the pressure on a solid is not 
necessarily hydrostatic. In any vanishingly small volume in a solid in which 
the pressure is non-hydrostatic, the pressure acting upon a plane passing 
through that volume depends upon the orientation of the plane. In the termin- 
ology of Harker (1939) the solid is under a “shearing-stress.” If non-hydro- 
static stresses are involved, the work done in any process is related not only 
to changes in volume (dilatation), but also to changes in shape (distortion). 
The pressure-volume terms in equations (1) and (2) must be replaced, in 
this event, by a more complicated expression taking into account the vectorial 
properties of the stresses and the resulting strains. For simplicity, however, 
we may consider the free energy of a substance under non-hydrostatic stress 
by comparing it with that of the same substance, at the same temperature and 
mean hydrostatic pressure, but not under non-hydrostatic stress. The mean 
hydrostatic pressure, P,,,, is defined as the arithmetic mean of the three princi- 
pal stresses. We may then, for convenience, define a quantity, D, by the 
following equation 
D = Gu-G, 
(15) 


where Gy is the free energy of the stressed solid and G, is the free energy 
of the same solid, unstressed, but at the same temperature and acted upon by 
a hydrostatic pressure equal to the mean hydrostatic pressure acting upon 
the stressed solid, 

A theoretical investigation of the probable magnitude of D, compared 
with the other factors affecting the free energy of a solid, has been made by 
Verhoogen (1951) who concluded that in general it should be negligible. To 
the extent that Verhoogen’s conclusion is correct this means that equations 
(1) and (2) may be applied to solids under non-hydrostatic stress as well 
as to substances under simple hydrostatic pressure, if “P” represents the 
arithmetic mean of the three principal stresses, when the pressure is non- 
hydrostatic. This means, further, that at a given temperature and mean hydro- 
static pressure, the free energy of a solid should be independent of “shearing 
stress,” a conclusion obviously in conflict with Harker’s concept (1939) of 
“stress” and “antistress” minerals. This concept implies that D is not only 
significant but varies significantly from mineral to mineral, being relatively 
large (if positive) in antistress minerals and relatively small (if positive) in 
stress minerals. Experimental evidence (Larsen and Bridgman, 1938) on 
the matter is sparse and negative, in that no significant effect of non-hydro- 
static stress upon chemical equilibria has yet been demonstrated, other than 
what can be related to the effect of the mean hydrostatic pressure involved. 

A few further comments might be made. One is that it is to be anticipated 
that D, if other than zero, should, when the crystal is anisotropic, depend 
upon the orientation of the crystal relative to the principal stresses. Such 
variations, furthermore, might be fully as great as those between different 
minerals, as, for example, are variations in hardness and optical properties. 
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A consequence of this is that stress-sensitive anisotropic minerals should 
show preferred orientation when formed under non-hydrostatic stress. Actually 
the most commonly cited anisotropic stress (kyanite, staurolite, chloritoid) or 
antistress (andalusite, cordierite) minerals show preferred orientation less 
commonly than do minerals such as mica, chlorite and quartz which are gen- 
erally regarded as relatively insensitive to non-hydrostatic stress. 

A second point is that the occurrence of a mineral in a highly deformed 
rock does not necessarily mean that it grew under strong non-hydrostatic 
stress, even where it can be demonstrated conclusively that the mineral formed 
during the deformation. Extreme deformation might as well be interpreted as 
evidence of extreme plasticity as of a strong non-hydrostatic stress. The more 
plastic a material is, in fact, the less capable it is of supporting a non-hydro- 
static stress. It is also true that strong non-hydrostatic stresses may be set up 
without large scale deformation, as, for example, when a crystalline aggregate 
is subjected to a sudden change in temperature. This might happen in thermal 
metamorphism at an igneous contact. Such stresses would, of course, differ 
in orientation from grain to grain but might nevertheless be great enough to 
bring about rupture of relict minerals, and are supposedly responsible 
(Harker, 1939) for the decussate (criss-cross) texture of hornfelses. In short, 
it would seem that the assumptions that the minerals of thermally metamor- 
phosed rocks grew in the absence of non-hydrostatic stress, and that those in 
the regionally meiamorphosed rocks of orogenic belts grew under relatively 
strong non-hydrostatic stress, are open to question. All considered, the writer 
is inclined to believe that the burden of the proof lies with the stress-mineral 
concept. The catalytic effects of non-hydrostatic stress are, of course, well 
known and have been discussed at some length by Eskola (1934). Significant 
displacements of chemical equilibria, however, are still unproven, and the 
writer proposes to show in the following discussion that many of the effects 
commonly attributed to non-hydrostatic stress are amenable to alternative 
explanation. 

In accord with the above remarks, then, we shall return to equations (1) 
and (2), bearing in mind, owing to our neglect of non-hydrostatic stress, 
that under some conditions they may not be entirely adequate. 

Equilibria involving a fluid phase—If a fluid phase is involved in the 
reaction we cannot assume constancy of AS and AV unless the fluid is known 
to be a liquid or in a liquid-like super-critical state. The thermal expansion 
and compressibility of most liquids is small and their heat capacities do not 
differ greatly from those of the corresponding solids. Slopes of equilibrium 
curves on a P-T diagram are likely to be steep (at a high angle to the tempera- 
ture axis) relative to those of solid-solid transitions owing to the considerable 
increase in entropy in passing from the solid to the liquid state. The volume 
changes per gram-atom are comparable to those in solid-solid transitions. The 
slopes are generally positive although the ice-water transition is an obvious 
exception. 

If a gas is involved, however, thermal expansion and compressibility are 
no longer negligible and AC, may also be large. Owing chiefly to the large 
molar volumes and high compressibilities of gases at low pressures, an 
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equilibrium curve involving production of a gas phase generally has a gentle 
positive slope at low pressures. With rising pressure the curve becomes both 
steeper and straighter as the gas becomes more like a liquid in its properties. 
Several equilibria of this type have been investigated experimentally (Bowen 
and Tuttle, 1949; Yoder, 1952). The slopes at the highest pressures at which 
investigations have been carried out (2-3000 bars) are steep, relative to most 
solid-solid transitions, and apparently positive (table 4), but the data are not 
sufficiently detailed to determine the slopes in the high-pressure region within 
very narrow limits. 


TABLE 4 
Slopes of Dehydration Curves at High Pressures 


Reaction Slope (bars/deg) 








Brucite + Serpentine —> Forsterite + Vapor 53 to 210 
Serpentine —» Talc + Forsterite + Vaper 49 to 240 
Tale + Forsterite —» Enstatite -+ Vapor 26 to 103 
Tale — Enstatite + Quartz + Vapor 21 to 68 
Brucite — Periclase + Vapor —34 to 34 
Chlorite — Forsterite + Cordierite + Spinel + Vapor 26 to 52 








Slopes are from experimental data of Bowen and Tuttle (1949) and Yoder (1952). The 
range represents the greatest and least possible slopes consistent with runs at pressures 
greater than 5000 PSI (350 bars). 


There is, however, reason to question the applicability of these curves to 


analogous reactions taking place in an actual geologic environment. We cannot 
assume that naturally occurring fluid phases, when present, are “pure com- 
ponents,” nor can we assume that they were subjected to the same mean 
hydrostatic pressure as that acting on the adjacent solids, Pure or nearly 
pure water is a possibility, but pure carbon dioxide is less likely, and we 
cannot assume that the molar properties of pure carbon dioxide in any way 
approximate the partial molar properties of carbon dioxide in a dilute aqueous 
solution at the same temperature and pressure. The chemical potential (partial 
molar free energy) of any component in a solution must be less than the molar 
free energy of the pure component at the same P and T. Otherwise, the com- 
ponent in question would separate as a distinct phase, the system thereby 
attaining a lower free energy. If a stable solution exists it is thus evident that 
the chemical potentials of the various components are lower than the molar 
free energies of these components were they to occur as pure phases at the 
same pressure and temperature. This means that if the “volatile” component 
does not occur as a pure phase, the actual transition temperature at any given 
pressure will be lower than that given by the experimental results mentioned 
above. It can be shown that, as the mole fraction or concentration of a given 
component approaches zero, its chemical potential must decrease indefinitely. 
As the concentration of the critical component in the fluid approaches zero, in 
other words, the equilibrium temperature must approach absolute zero, what- 
ever the pressure. 

Even if the critical component does occur as a pure fluid phase its molar 
free energy will decrease, at any given temperature, as the pressure upon it 
is lowered. This means that even where the critical component can be assumed 
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“pure” we will overestimate the transition temperature from the above curves 
if the pressure upon it is less than that on the adjacent solids. 


It is evident, then, that the curves discussed above are useful only as 
limiting cases, and tell us for a given “rock-pressure” only the maximum 
possible transition temperature. They are applicable only where the fluid con- 
tains only the “volatile” component in question and is acted upon by the same 
pressure as the adjacent solids. 


A further complication arises in that we cannot assume that a true fluid 
phase is always present when reactions of this type take place in actual rocks. 
Geologic and petrographic evidence suggests that in many instances the 
“volatile” or mobile components have been removed by diffusion as rapidly 
as produced, or supplied by diffusion as rapidly as consumed. True fluids 
of course exist in fissures and other large openings, but “volatilizations” and 
“devolatilizations” have clearly not been restricted to those portions of the 
rock in direct contact with such openings. The “intergranular film,” further- 
more, is believed to be quite unlike a liquid or gas in its properties. The water 
in any subcapillary opening is strongly bonded by adsorption to the adjacent 
surfaces and probably best regarded as solid — analogous, structurally, 
to zeolitic water or water of hydration. The intergranular film is also quanti- 
tatively inadequate to account for the vast quantities of water given off in 
the metamorphism of a shale to a sillimanite gneiss, or for the carbon dioxide 
produced in the conversion of a siliceous dolomite to a calc-silicate granulite. 
It appears, then, that a rock must be regarded as an open system, permeable 
to certain relatively mobile components such as water and carbon dioxide. 


OPEN SYSTEMS 
Partial molar properties.—Partial molar properties may be defined by 
equations of the type: 


_ aX, 
cae (5) parmeme one (16) 
where X may be any extensive property of phase a such as V, S, E or G, i 
refers to the component in question, and j to all other components of a. Thus, 
Xia is to be read: “partial molar X of i in a.” 

The partial molar Gibbs free energy, however, is commonly referred to 
as the chemical potential and denoted by the symbol », so that 





=) 
— (= P,Tnja,Mxa; eve (16a) 

For two phases sharing a given component i to be in equilibrium, »; in 
each must be the same. Otherwise the free energy of the system may be 
lowered by appropriate transfer of i. 

Significance of chemical potential in open systems.—Similarly if a 
system in thermal equilibrium with its surroundings is open with regard to a 
mobile component /, then »; must be the same in both the system and its 
immediate surroundings in order that there be no tendency for transfer of 
i to or from the system. The free energy of the “megasystem” consisting of 
the open system plus its environment is at a minimum with regard to distribu- 
tion of component i when such a distribution has been attained. It is interest- 
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ing to note that the megasystem as defined above is a closed system even 
though both its component parts (open system and environment) are not. 
The distinction between a closed and open system is thus dependent in large 
part on what one chooses to define as the limits of the system. Although closed 
systems are easier to treat thermodynamically, it is impractical in many 
instances to hunt about for a “megasystem” large enough to remain “closed” 
through a given process. We might be forced in certain instances to consider 
not only a large portion of the Earth’s crust but the hydrosphere or atmo- 
sphere as well. A system of such size is not only cumbersome but also rarely, 
if ever, in a state of complete internal equilibrium. For these and other reasons 
we shall find it simpler to regard rocks as open system with regard to certain 
mobile or “volatile” components. 

Multivariant equilibrium.—tlf a rock with a given content of fixed com- 
ponents is an open system with regard to a mobile component, 1, its most 
stable state (mineral facies) will be dependent not only on pressure and tem- 
perature but also upon the value of ; in its immediate environment. A state 
of equilibrium between two different facies may now be expressed by the 
equation 

AG =0= AEs T PAVs _ TASs + py An; (17) 
Where AG refers to the megasystem, rock plus environment, AEs, AVs and 
ASs refer to changes in the rock, and An, is the change in the number of moles 
of i in the environment (= — An;,). 

It can be shown that the chemical potential may also be defined by the 
relation 


_ fe 
Lia — — , 
Oniaf S,VnjasMxa; . - - (16b) 
so that the combined first and second laws for open systems have the form 


dE, = TdSs — PdVg + Sipidni, 


(18) 

and the differential form of (17) is then 
d(AG) = 0 = AVgdP — ASgdT + Anjdp; (19) 
Equations (17) and (19) are analogous to equations (1) and (2) for 
completely closed systems. Equation (17) is that of a surface of divariant 
equilibrium in a three-dimensional P-T-»; diagram. The field of stability of a 
specific mineral facies in such a diagram is a volume bounded by divariant 
surfaces (rather than an area on a P-T diagram bounded by univariant 
lines). The facies-volume is a region of trivariant equilibrium, and the maxi- 
mum number of phases in a facies is equal to the number of fixed components. 
The intersection of two divariant surfaces is a univariant line, and of three 
an invariant point. At the invariant point the number of phases may exceed 

by two the total number of components. 

Equations (17) and (19) may be generalized to 
AG =0=> AEg + PAVs ~— TAS; + Yip; An, (17a) 
and 

d(AaG) = 0= AVsdP - ASsdT + y,Anjdy; (19a) 
if two or more mobile components are involved, and they are now equations 
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of trivariant or multivariant equilibrium. The equations are no longer amen- 
able to simple graphical representation but we may think of the field of 
stability of a mineral facies as a region in (m + 2)-dimensional space limited 
by equations of (m + 1)-variant equilibrium, where m is the number of 
mobile components. Fortunately, few specific reactions involve more than 
one mobile component. For a random choice of P, T and all p; we are, in 
general, in a facies-region and the variance is m + 2. This means that in 
open systems of this sort Goldschmidt’s “mineralogical” phase rule should be 
modified to 
p=c-m 
(20) 


where p is the maximum number of phases, c is the total number of com- 
ponents, and m, as above, the number of mobile components. This relationship 
was first derived by Korzhinsky (1936, 1950). It has, moreover, commonly 
been assumed or implied in triangular diagrams for mineral assemblages, in 
that the water or carbon dioxide content of a phase is generally neglected. 
The fact that this procedure is successful has been taken in some instances as 
indicating that water or carbon dioxide was present in excess (Bowen, 1940; 
Yoder, 1952) during the metamorphic process. In view of the lack of other 
geologic or petrographic evidence that water or carbon dioxide was present as 
a separate phase in such instances, it appears more reasonable to draw the 
alternative conclusion, specifically, that a petrologic system is open to water 
and carbon dioxide. 

It will be noted that a system in which all components are mobile will, 
for a random choice of P, T and all »;, have “no” phases. This appears anoma- 
lous but is not, since we have little interest in systems without phases. Better 
ways of stating the situation might be to say that for a random choice of 
P, T and all »; there is, in general, no actual phase that will have the specified 
properties, or to say that the variables are not all independent. When a phase 
exists in which all components are mobile, the equilibrium between this phase 
and its environment is univariant, that is, for any random choice of P and 
all »;, T is specified, or for a choice of T and all p;, P is specified. Clearly no 
more than two phases containing only mobile components can ever coexist 
stably in an open system and these only at a specific temperature and pressure. 

The limiting cases.—The problem, then, is to consider the ways in which 
the chemical potential of a given mobile component, i, may vary within the 
crust of the Earth. Two extreme limiting cases have already been mentioned. 
To repeat, these occur: (1) when the chemical potential of i is allowed to 
decrease indefinitely (u; — — 0) and (2) when the chemical potential of 
i is allowed to approach its maximum possible value for a given pressure and 
temperature (#; — »imax.e.>). 

Regarding i as a component in a solution, condition (1) above is fulfilled 
as the mole-fraction (N,) of i in the solution approaches zero, or, for our 
purposes, when i is virtually absent from the immediate environment of the 
system in question. In this event no phase containing i is stable at any pressure 
or temperature. 
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Condition (2) is effectively fulfilled when the stable form of pure i is 
present (N; = 1), making pjmax. equal to the free energy per mole of the most 
stable form of pure i at the given pressure and temperature. Values of ; in 
excess of this can never represent the most stable state inasmuch as there must 
then be a decrease in free energy upon separation of the “pure i.” Condition 
(2) thus corresponds, in essence, to the condition of “excess water” (= ex- 
cess i) discussed by Yoder (1952). As long as condition (2) is fulfilled the 
equilibrium is univariant and essentially identical to that in a closed system 
containing excess i as a nearly pure phase. The qualification “essentially 
identical” enters in that condition (2) is defined with reference to “strictly” 
pure i, whereas actual closed systems contain at best “nearly” pure i. This 
means that conditions (1) and (2) also represent the effective limiting values 
of »; in closed systems. In a closed system, however, it is nj, not »; that must 
be specified to determine the state of the system, and the two variables are by 
no means equivalent. Specifically, the amount of a phase containing i may be 
changed at P, T without affecting »;, but obviously changing nj. Similarly, 
the statement that a hydrate and its anhydride are in equilibrium at P, T fixes 
enzo, though nyxo may be changed by simply varying the amount of either 
phase, and conversely if component j is removed at P, T from a solution con- 
taining i, »; will be changed but n, will not. 

Distribution of the mobile components.—There are at least two alterna- 
tive hypotheses which might be considered to represent the stable vertical 
distribution of a given mobile component, i, in the crust of the Earth. One is 
to assume that the distribution of i is such that there is no longer any tendency 
for vertical diffusion of i. Complexities arise here, however, owing to the 
difficulty in determining the diffusion potential in non-isothermal systems. 

The other hypothesis, which we shall investigate first, is to suppose that 
the rock is in osmotic equilibrium with an adjacent fissure system containing 
a dilute aqueous (hydrothermal) solution at the same temperature. Here, we 
shall consider the closed megasystem (rock plus nearby fluid) rather than 
the rock alone. It is instructive, however, to review first the general principles 
of osmotic equilibrium, 

Osmotic equilibrium.—Perhaps the simplest example of osmotic equili- 
brium is that in which the system is divided into two parts separated by a 
partition or membrane, permeable with regard to certain components, imper- 
meable to others. When equilibrium is reached under such conditions, it is 
found experimentally that the hydrostatic pressure upon the two parts of the 
system is different, the difference being supported by the semi-permeable 
partition. Conversely, the distribution of the mobile constituents relative to 
the partition is found to depend on three independent variables: the tempera- 
ture (which at equilibrium must be constant throughout the system), and the 
pressures on each side of the partition. As an example let us consider an 
equilibrium of the sort: 

hydrous solid = anhydrous solid + water 
and let us suppose that the solid is separated from the fluid by a wall per- 
meable to water, and that the pressures-on opposite sides of the wall may be 
different. We have then, at equilibrium 
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If the “fluid” is not pure water such a transition will take place over 
an interval owing to the changing composition of the fluid. If the fluid is 
kept at constant composition, possibly through a continuous supply of “fresh” 
fluid, however, the transition will be abrupt, and if the volume of fluid is 
large relative to the reacting solids any interval must be narrow, since the 
mole fraction of water in the fluid cannot then change appreciably in the 
course of the reaction. In any event we may consider, alternatively, a meta- 
stable equilibrium within the actual interval, as outlined in an earlier section 
of this paper. 

From equation (1) we may also write 

AG = 0 = AE-— TAS + PsAVs + PpAVe (22) 
and from equation (2) 
d(AE) = Td(AS) — Psd(AVs) — Prd(AVr) 
(23) 


hence, if equilibrium is to be maintained, 
d(AG) =0=> AVsdPs + AVrdPr — ASdT (24) 
It must be emphasized here that Py is the actual pressure on the fluid, 
and not to be identified with its vapor pressure, if liquid, nor with the “par- 
tial pressure” or “partial vapor pressure” of any of its components. Holding 
each of our variables Pr, Ps and T constant in turn, and maintaining equili- 


(i): -2 
oT Pp AVs 


brium we have 


(25) 


(2h 
oT P. AVy 


(3) 
aPs 7 wee) 


A useful relation? is 


( 7) AS 
aT esa 

G=0 AV,+AVr —~ 

A 8 * dPs (28) 
The relative fields of stability of rocks A and B in the above example can be 
represented by a three-dimensional diagram based upon the variables Pr, 
Ps and T, the surface separating their respective fields of stability being one 
of divariant equilibrium. It is instructive, however, to consider first equations 
(25), (26) and (27) above. 


Equation (25) is the slope of the trace of the equilibrium surface on a 
section at constant Py. AS is defined as positive and AVs, with few exceptions 
(table 5), will be negative, hence it may generally be assumed that the curve 
will have a negative slope. Its curvature will depend chiefly upon variations 


? With the stipulation Py = Ps = P equations (22) and (28) become identical with 
equations (1) and (3) for ordinary, non-osmotic equilibrium. 


(26) 
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in AS, but will be slight if ACp is small, that is, if the fluid is a liquid or in 
a liquid-like supercritical state. 

Equation (26) is the slope of the corresponding curve on a section at 
constant Ps. AVp can be negative only if the mobile component in question 
has a negative partial molar volume in the fluid. This is possible but unlikely. 
This curve should thus have in general a positive slope, and, particularly if 
the fluid is a gas, be markedly concave toward the Py-axis. The curvature will 
be most pronounced at low Pp where the compressibility of the fluid is high. 

Equation (27) is the slope of the corresponding curve on a section at 
constant temperature. It should, in general, have a positive slope, and if the 
fluid is a gas, be markedly concave to the Pg-axis in regions of low Py. 

Another section of some interest is that defined by the stipulation that 
Py = Ps. The curve on this section is essentially that for an ordinary non- 
osmotic equilibrium and will have the form of the experimental curves dis- 


cussed earlier, though offset toward lower temperatures if the fluid is not 
“ ” 
pure. 


TABLE 5 
Change of Volume with Dehydration 


Reaction AVs (cc) (cc) 


Avs 
Anueo 
Brucite —» Periclase + Water -13.0 -13.0 
Serpentine — Forsterite + 2 Water + Enstatite -33.5 -16.8 
Tale — 3 Enstatite + Quartz + Water -19.1 11.1 
Phlogopite + 3 Quartz —» 3 Enstatite + Orthoclase + 

Water -17.7 -17.7 
2 Diaspore —>» Corundum + Water -9.2 -9,2 
2 Gibbsite —-» Corundum + 3 Water —38 -12.7 
Kaolinite —» Kyanite + Quartz + 2 Water -32.6 -16.3 
Pyrophyllite —» Kyanite + 3 Quartz + Water -16.8 -16.8 
Pyrophyllite —» Andalusite + 3 Quartz + Water -8.8 -8.8 
Pyrophyllite —> Sillimanite + 3 Quartz + Water -118 -11.8 
Muscovite + Quartz —» Kyanite + Orthoclase + Water -15.4 
Muscovite + Quartz — Andalusite + Orthoclase + Water -7.4 -7.4 
Muscovite + Quartz —> Sillimanite + Orthoclase + Water -10.4 
Chlorite + 2 Quartz —» Pyrope + 2 Enstatite + 4 Water -20.0 
2 Chlorite + 7 Quartz —» Cordierite + 8 Enstatite + 

8 Water 
3 Chlorite + Muscovite + 3 Quartz —» 4 Pyrope + 

Phlogopite + 12 Water -20.0 
3 Staurolite + 2 Quartz —» Almandine + 5 Kyanite + 

3 Water -14.3 
3 Staurolite + 2 Quartz —» Almandine + 5 Andalusite 

+ 3 Water : -1.0 
3 Staurolite + 2 Quartz —» Almandine + 5 Sillimanite 

+ 3 Water -18 -6.0 
Analcite — Jadeite + Water -38 ~38.0 
2 Analcite —» Nepheline + Albite + 2 Water -37.5 -18.8 
Lawsonite —> Anorthite + 2 Water -3 -1.5 
4 Clinozoisite + Quartz —» Grossularite + 5 Anorthite 

+ 2 Water +44.6 +22.3 
2 Clinozoisite + Muscovite + 2 Quartz —» 4 Anorthite 

+ Orthoclase + 2 Water + 62.6 +31.3 
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Fig. 1. Probable form of the surface of divariant osmotic equilibrium for reactions 
of the type: Hydrate + Fluid -» Anhydride + Fluid. 


The probable appearance of such an equilibrium surface is shown in 
figure 1. Clearly we have, here, no interest in the region where Py exceeds Ps 
and may take as our limiting cases the conditions: Py = 0 and Py = Ps.* 

The special case of equation (25) where Py is constant has been given 
by Johnston and Niggli (1913) and later by Harker (1939) as representing 
the effect of increasing pressure upon a transition temperature when that 
pressure increase acts only upon the solid phase. For equation (25) to be 
applicable to equilibria in the Earth’s crust, however, a situation must exist 
wherein the rock pressure is variable but the fluid pressure remains constant. 
Inasmuch as both rock and fluid pressure may be expected to vary there does 
not seem to be any justification for the use of equation (25). 


EQUILIBRIUM BETWEEN A FLUID-FILLED FISSURE 
AND THE SURROUNDING ROCK 


A situation in which the principles of osmotic equilibrium might legiti- 
mately be applied is in the chemical equilibrium between a column of rock 
and a column of fluid in a nearby fissure. Obviously such a model cannot 
represent a state of hydrostatic equilibrium unless rock and fluid are of the 
same density, which is highly unlikely. The immediate walls of the fissure 
in such a model must be under a non-hydrostatic stress, and must also be 


* The use of Py rather than yw; as the extra variable in dealing with osmotic equilibrium 
is chiefly a matter of convenience, since Py and y; are not independent in a fluid of 
given temperature and composition. The condition Py = Ps corresponds, where the fluid 
is pure i, to pimex.. If Py were to exceed Ps the free energy of the system would be 
lowered by the passage of i through the membrane so as to form a phase of composition 
i on the “solid” side of the membrane. Py greater than Ps thus does not correspond to an 
equilibrium state where the membrane is permeable to all fluid components. 
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subjected to a mean hydrostatic pressure considerably less than that prevail- 
ing in the rock as a whole at the same depth, Because of this, the chemical 
equilibrium between the fluid and the portions of the rock immediately ad- 
jacent to it will be different from that between the fluid and those portions 
of the rock at sufficient distance from the fissure that the mean hydrostatic 
pressure may be regarded as determined solely by the weight of the overlying 
column of rock. For the moment, then, we will set aside the problem i the 
chemical equilibria within the fissure walls and regard the “wallrock” simply 
as the semi-permeable partition supporting the difference in pressure between 
rock and fluid. 

Hydrostatic equilibrium.—Before proceeding with the fissure system we 
must first consider the vertical variation of pressure in any medium at hydro- 
static equilibrium in a gravitational field. Under these conditions 

dP = pgdz (29) 


where p is the density, g the acceleration of gravity and z the depth taken 
as positive downward. At any point in such a column 

P-P*= f,"pgdz (30) 
where P* is the pressure when z = a. Equation (30) can also be written 


P—P* = p"g (z-a) (30a) 


where p™ is the mean density of the column between z and a. It is generally 
safe to assume that mean rock density, p,", is about 2.7 gms/cc. 

For the fluid column it is difficult to predict the value of pp” between 
any two points in the column or even the value of pp at any given point. In 
general, however, we are perhaps safe in assuming that the fluid is a dilute 
aqueous solution, closely resembling, in its physical properties, pure water 
at the same pressure and temperature. Fortunately, the density of pure water 
is known over a wide range of temperatures and pressures (Kennedy, 1950). 

It is convenient to make use of the relations 

_ Mr 
Pr = 


Vr (31) 


(where My and Vy are the mean molecular weight and mean molar volume, 
respectively, of the fluid at the point in question) and 


(32) 
where the superscript signifies the corresponding mean quantities over the 
interval from a to z. 

We now have: 
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Eliminating the terms containing z or dz we may now express Py and 
dP» in terms of Ps and dPg so that 
M 
dPy = ™ dPs =—= Ps 
Ps psVr 


or, since local variations in rock density are of no concern here 
Mp 
dPy =——= dPs 


Ps F 


Py — P;* — PP (Ps _ Ps") — A (Ps - Ps*) 
Ps ps" Ve" 
(37) 
If our fluid column extends to the Earth’s surface it is convenient to 
have a at zero depth so that Pp* = Ps* = atmospheric pressure (1.013 bars). 
In the event that the fluid column does not extend to the surface, or is ob- 
structed by some relatively impermeable barrier, it is convenient to have 
a represent the top of the free fluid column. In this case Py“ is equal to or less 
than P;*. 
With the assumption that our fluid is a dilute aqueous solution, we may 
substitute the molecular weight of water (18 gms) for My and Mp” so that 
our final equations of hydrostatic equilibrium become 


ve 
F Vr s 


Py _ P;* = " (Ps ~ Ps") 


where V° is a constant with the dimensions of volume defined by 


v= ~ = =. = 6.7 ce. 
Ps Ps 

(40) 

For a column open to the Earth’s surface, and at a depth such that Pp* and 


P;* are negligible compared with Pp or Ps, we may use the relation 
* 


Py == Ps 


Vr 
in place of equation (39). (41) 
The quantity leading to the greatest uncertainty in the above relations 
is Vp", the mean molar volume of the overlying fluid column. For some 
purposes Vp" may be assumed constant at about 18 cc/mole, the molar 
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volume of water at 25°C. and atmospheric pressure. In this case the ratio 

* 
— is = and equations (39) and (41) may be simplified accordingly. 
If the temperatures and pressures over a significant portion of the fluid column 
are such that this assumption is unreasonable, Ve" may be estimated from 
the tables of Kennedy (1950). 

Osmotic equilibrium in a fissure system.—The equations of osmotic 

equilibrium for the model fissure system now take the form: 

y* y* 
. Tv." + Ps ZF aAVe . ras (42) 


dP i as 
(sr) AG =0 “i 2 v* 
AVs + = AVy 
Vr 


AG = 0 = AE + PsAVs + (Pr*—P 
and 





(43) 
The practical use of these equations is hampered by the difficulty in 
evaluating the quantities AVp. Ve". AS and AE. If, however, we restrict our 
discussion, at first, to reactions where the only volatile component involved 
is the solute itself—in this case, water—we may make simplifications not 
otherwise possible. It is fortunate in this regard that many of the petro- 
logically important reactions are simple hydrations or dehydrations. 
Reactions with water the only mobile component.—In a sufficiently dilute 
solution the partial molar properties of the solute, the corresponding mean 
molar properties of the solution as a whole, and the corresponding molar 
properties of the pure solute, at the same pressure and temperature, are vir- 
tually identical quantities. If the property in question is volume this may be 
expressed formally by the equation 


Vive = Vr = Vir‘ (P,T) (44) 


where the subscript, i, refers to the solute and the quantities are given in the 
order mentioned above. 

Relations analogous to (44) also hold for other extensive properties 
such as S, E, H, and G except that conventional notation in the case of the 
Gibbs free energy leads to 


hir = Gr - Pir (45) 


Because of (44) we then have, where i is the solute, 


AVy — VrAnir (46) 


and equation (43) may immediately be simplified to 





dP = AS 
(5 AG 0 ais AVs oo V’Anip (47) 


To simplify equation (42), however, a further assumption is necessary, 
specifically 


(48) 
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that the mean molar volume of the fluid at the point in question is nearly 
equal to its mean value taken over the entire overlying column. This is cer- 
tainly valid where the column is in a liquid or liquid-like, supercritical state 
throughout, a condition probably fulfilled in most cases. The range of validity 
of this approximation, however, is actually somewhat greater than this be- 
cause the increase of pressure with depth tends to offset the effects of the 
rise in temperature, and (48) is thus valid wherever the thermal gradient 


oT 
in the column approaches (Sp ) V for water. With approximation (48) 
equation (42) simplifies to il 
AG = 0=> AE+ PsAVs + (Pp"*Ve —Ps‘V P +PsV*) Anjp — TAS (49) 


From equation (41) (valid where P*s and P*y are negligible), this in turn 
becomes 

AG = 0 = AE + Ps(AVg + V*Anig) — TAS (50) 
The close analogy between equations (47) and (50) and equations (3) and 
(1) respectively is interesting, V° is evidently to be interpreted as the effective 
molar volume of water under the above conditions of “fissure equilibrium.” 
It is also instructive to compare equation (47) with equation (25). 

From the values in table 5 it is evident that AVs per mole of water lost 
is generally on the order of 8 to 20 cc/mole of water, The values for reactions 
involving epidote or lawsonite versus plagioclase are exceptional, probably 
because the dehydration is incidental, here, to major crystal-structural changes 
analogous to those in the reaction 2 jadeite > nepheline + albite. 

Entropy determinations on hydrated minerals are few, the only data 
currently available being for the hydroxides of magnesium, aluminum and 
calcium. The reaction 2 gibbsite — corundum + 3 water involves an entropy 
change, at 25°C. and one atmosphere, of 408 + 18 decijoules/deg/mole of 
water, and the reaction brucite — periclase + water an entropy change of 
344 + 21 decijoules/deg/mole of water. Noting that these dehydrations do 
not involve a change in the first-order coordination of the aluminum or 
magnesium atoms, it might be suggested that the mean of these values (375 
dj/deg/mole of water) is a fair approximation wherever this condition is 
fulfilled, and the water is held as hydroxyl. This assumption, though ad- 
mittedly shaky, has been employed in compiling table 6, for reactions in which 
there is no major change in the first-order coordination of the metallic ele- 
ments, Where there is a coordination change, as in the breakdown of pyro- 
phyllite to andalusite or sillimanite, rather than to kyanite; or in the break- 
down of analcite to jadeite, rather than to nepheline and albite, the AS has 
been modified in accord with the data in table 3. It will be observed that the 
values of dP/dT so calculated are negative and somewhat steeper than the 
curves for most solid-solid univariant equilibria. In any region where the 
fluid is a gas or gas-like supercritical fluid AS will certainly be greater than 
our estimate and the slopes of the curves steeper than in table 6, If the mean 
density of the overlying fluid column is significantly less than that at the 
point in question, moreover, Vy" and the dehydration temperature at any Ps 
will be less than predicted by relation (48). Dehydration temperatures under 
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the conditions of “fissure equilibrium” should thus, in general, vary with depth 
(Ps) as indicated by curve ¢ in figure 2. Curve a is the general form of curves, 
calculated or determined experimentally, where Py is equal to Ps and the fluid 
is pure water. Curve 6 is that predicted by the Johnston-Niggli relation (25), 
and curves d and e suppose “fissure equilibrium” beneath an impermeable 
barrier. 

Metamorphism of carbonate rocks.—It must be emphasized that the 
above simplifications cannot safely be made in dealing with reactions in- 
volving carbon dioxide. Specifically, we cannot make any simple assumption 
as to the partial molar volume, or partial molar entropy of carbon dioxide 
in the fluid, even with the assumption that it is a dilute aqueous solution. The 
mixing properties of water and carbon dioxide, owing to the formation of 
carbonate and bicarbonate ions, are complex and are sensitive to the presence 
of other components in solution. It is clear that curves such as that calculated 
by Goldschmidt (1912) represent maximum decarbonatization temperatures 
at any depth, but we cannot, at present, evaluate equations (42) and (43) 
with sufficient confidence to say how much the conditions of fissure equili- 
brium would lower these temperatures in any specific instance. 

Goldschmidt’s curve for the calcite-wollastonite equilibrium has been 
recalculated recently by Danielsson (1950). Danielsson also, in the same 
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Fig. 2. Probable forms of dehydration curves assuming the fluid to be a dilute 
aqueous solution: 

(a) Py = Ps 

(b) Py = constant (Johnston-Niggli equation) 

(c) The fluid occupies a fissure system open to the surface 

(d) The fluid occupies a fissure system beneath an impermeable barrier at the 
base of which Py = Ps 

(e) The fluid occupies a fissure system beneath an impermeable (or semi- per- 
meable) barrier at the base of which Pr < Ps 

Curves (c), (d) and (e) are uniquely determined for a given dehydration only 

in the event that the mean density of the overlying fluid column is known. 
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paper, considered the reaction from the point of view of simple osmotic equi- 
librium, assuming the fluid to be pure carbon dioxide. He pointed out, how- 
ever (p. 69), that the presence of water vapor and, implicitly, other com- 
ponents would alter the results. 

The walls of the fissure—The equilibria within the fissure walls, as 
distinct from those in the main body of the rock, constitute an interesting 
problem in view of the extensive current literature on “wallrock alteration.” 
The stress in the walls must be non-hydrostatic wherever Py is less than Ps, 
the mean hydrostatic pressure on the solids, P,", tending to approach Py 
near the fissure. Assuming that the unequal stress distribution affects the 
equilibrium only by way of P,”, it is probable from the equations of osmotic 
equilibrium that we should anticipate higher dehydration temperatures in 
the immediate vicinity of the fissure than at a distance. In other words we 
may expect the degree of hydration of a given rock to increase as the fissure 
is approached, other things being equal. Conversely, we are not justified in 
concluding from an observed phenomenon of this sort that “other things” 
were not equal. It does not, for example, necessarily mean that the fluid in 
the fissure was at a lower temperature at the time of alteration, inasmuch 
as the formation of a fissure filled with fluid could affect the wallrock in 
this way even if the fluid were slightly hotter. If the degree of hydration de- 
creases toward the fissure, vein or dike, however, it is probably safe to assume 
that the material in the fissure was significantly hotter than the surrounding 
rock. In short a “contact metamorphism” gives us positive evidence of a 
thermal gradient whereas a “wallrock alteration” does not. 


VERTICAL DIFFUSION OF THE MOBILE COMPONENTS 

An alternative hypothesis is to suppose that an equilibrium distribution 
of the mobile components in the Earth’s crust can be effectively maintained 
by diffusion alone. In this view the chemical potential gradients by which we 
might define the limits of stability of a mineral facies would be those obtain- 
ing when any tendency for net transport of the mobile components by diffusion 
has ceased to exist. 

Isothermal diffusion—In a system at constant temperature throughout, 
net horizontal transport of any mobile component, i, will cease when the 
horizontal components of its chemical potential gradient are zero. That is, if 
the x- and y-coordinates are in a horizontal plane, 


Op Oni 
i, -\8),0° 
as die Be (51) 


when there is no net horizontal flow. Vertical flow will cease, however, only 
when the vertical component of the chemical potential gradient is balanced 
by the force exerted by gravity upon one mole of i, hence, at equilibrium 


Opi ane 
“Oz — Mig 
1 (52) 


where z is the depth, positive downward, M, is the mass of one mole of i, and 
g the acceleration of gravity. 
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Making use of equation (29), then, and neglecting any horizontal com- 
ponents, we have at equilibrium: 


(2) —_ M 
dPs} > ps” (53) 


The usefulness of the above relations to petrology is impaired by the 
constant-temperature restriction. Diffusion at igneous contacts is obviously 
not isothermal, nor is diffusion over any extended vertical range, owing to 
the geothermal gradient.‘ 

Non-isothermal diffusion—No system in which there is a net flow of 
matter (diffusion) or of heat may be regarded as in a state of true thermo- 
dynamic equilibrium. Since no known substance is a perfect thermal insulator, 
there can be no true thermodynamic equilibrium in any system in which 
there is a thermal gradient. 

There has, however, been an extension of thermodynamic methods to 
the point where such problems may be handled at least qualitatively, provided 
such a flow of heat or matter has reached a steady state (does not change with 
time). Recent summaries of this work are those by de Groot (1951) and 
Denbigh (1951). These researches are pertinent in the present instance be- 
cause we are interested in just such a steady state,—that in which there is no 
flow of matter and a constant flow of heat. The flow of heat may be held 
constant by some external agency, the nature of which need not concern us. 
The phenomenon of a compositional or chemical potential gradient resulting 
from a steady thermal gradient is known as the Soret effect. 

Following the method outlined in Denbigh (1951, p. 14-21 and 63-75) 


we arrive at the relation 


* 
grad p; = F, - (5 +e grad T (54) 
for the condition that there be no flow of matter in a system with a steady 
thermal gradient. F; is any external force acting on i (in this case gravity) 
and Q;* is a quantity known as the molar “heat of transport” of component 
i. Q;* may be defined by the relation 


Q\* _ Q, sa H, (55) 


where H, is the partial molar enthalpy of i, and Q, is the energy transported, 
per mole of i, in isothermal diffusion at P, T. This suggests correlation of 
Q\* with the extra energy possessed by the moving atoms as compared to 
the “average” atoms, and suggests, further, a close relationship between Q;* 
and the activation energy of diffusion, Q;* will presumably vary with tempera- 
ture, pressure and composition, but there is virtually no quantitative informa- 


* What is roughly equation (52) above has been derived by Ramberg (1944) and has 
appeared more recently in the book by Barth (1952). The necessary stipulation of con- 
stant temperature, however, is missing in Ramberg’s relation, as is also the acceleration 
of gravity. This last omission is valid only if it be understood that z is the negative of 
the geopotential rather than the depth, or, alternatively, if the forces and energies are 
understood to be in gravitational rather than absolute units, Neither condition is con- 
ventional usage, nor was either specified by Ramberg. The constant temperature stipula- 
tion has also been neglected in the discussion by Turner and Verhoogen (1951, p. 40-42). 
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tion available. From the probable relationship to activation energy, however, 
it appears reasonable to predict that Q,* will be a positive quantity at least 
in the case of solid or liquid diffusion. 

Neglecting any horizontal components of the geothermal gradient, and 
assuming hydrostatic equilibrium in the crust, equation (54) becomes 


dy; =P i Q,* dT 

ss (s a \= (56) 
or 

dy, = gp -(s + @ )ar (56a) 

ps T 
The integrated form of equation (56a) is 
— s Mi s r = Qi* Gh 
fi =< pi +— (Ps — Ps ) = (S; + _— )dT (57) 
ps a" T 


where the superscript a refers to values at the surface of the Earth or at the 
base of any impermeable barrier. Where a is at the Earth’s surface Ps" may 
be neglected and we have 
ja ™ P 
a= M - j (S; + Q: )dT (57a) 
Ps vw -," T 
Multivariant equilibrium.—F or simplicity we may define V;* and S,* by 


the relations 
Opi M, 
\ ” _— ( = : 
oP T ps” 


Opi _ * 
so (H), = (s+ 
Tp T (59) 


for steady state conditions with no net flow of i. It is apparent that these 
quantities are the effective molar volume and entropy, respectively, of com- 
ponent i in the immediate environment of a given open system under the 
above conditions. 

By substituting for »; and dy; in equations (17a) and (19a) we have 
then 


AG = 0= AEs T PsAVs — TASsg + +, An;[pi* + V,* (Ps — Ps*) - 








(58) 


and 





SaaS, * dT} 
(60) 
and 
dPs ee 
(a). 0 AV, + 3, V;* An, (61) 


In the absence of impermeable barriers equation (60) becomes 
AG = = AEs + PsA Vs — TASs = x >, An; (p;" + V,*Ps- SaaS; *dT) 
(60a) 
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Where water is the only mobile component we have then 
4G =-0= AEs + PsAVs - TASs + Anyeo( "Heo + . 
V* w0Ps - JS Ta 5S H20dT ) 
(62) 


and 





(st) a ASs + S* x20 Anno 
aT AG = 0 AVs + V*i20 Anuso 
(63) 


These may be compared with equations (50) and (47) for fissure equilibrium. 


It is likely that the partial molar entropy of “interstitial water” (Suzo) 
would be significantly greater than the effective molar entropy of the com- 
bined water in a given reaction { —ASs owing to its less ordered state 

( Anr0 
of aggregation. It is also probable that Qy2o* is a positive quantity if we are 
correct in relating it to the activation energy for diffusion. These considera- 
tions suggest that equation (63), like equation (47), would indicate steep 
negative slopes for most dehydration curves, and consequently that “diffu- 
sion equilibrium” and “fissure equilibrium” do not yield greatly different 
results. 

Diffusion versus flow through fissures —The merit of “diffusion-equili- 
brium relative to “fissure equilibrium” depends to a considerable extent upon 
the nature of the geologic situation to which applied. Where a fissure system 
exists the water in the fissure would presumably exercise the dominant control 
over the chemical potential gradient of water, at least, in the immediately 
adjacent rocks. 

For other mobile components the situation is less clear. It is possible. 
if the fluid in the fissure were not in active circulation, that the chemical 
potential gradients of the various solute species might reach equilibrium 
within the fluid column. These gradients would probably differ somewhat 


from those in the adjacent rock column owing to the dependence of S; and 
Qi* upon the environment of i, even where rock and fluid are in thermal 
equilibrium with each other at all depths. In this event the gradients estab- 
lished in the fluid column, when such is present, would probably dominate 
owing to the faster diffusion in the fluid. 

An opening where Py is less than Pgs is thermodynamically metastable, 
if not unstable, and may be expected to close in time, either through filling 
by solution and recrystallization or by mechanical collapse of the walls. The 
continued existence of a fissure system at any depth thus depends on its being 
maintained by some outside agency such as diastrophic or orogenic move- 
ments, or explosive igneous activity. It would appear, then, that no definite 
statement can be made as to the depth at which a fissure system open to the 
surface can exist, the answer depending upon such imponderables as the rate 
of closure of a fissure, once formed, versus the rate at which it is reopened by 
tectonic or igneous activity. A fluid-filled opening not open to the surface can 
exist at any depth as long as Py is not too much less than Ps, owing to support 
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of the walls of the fluid. In general, then, we might expect the relative im- 
portance of vertical solid diffusion to be greater for “deep” rather than 
“shallow” processes, partly because of the relative scarcity of openings at 
“depth” and partly because of the greater effectiveness of diffusion at high 
temperatures. 


GEOLOGIC APPLICATIONS 

We have considered several limiting cases with regard to reactions in 
which mob:le components are directly involved. In each of these cases, we 
have, in effect, treated the chemical potential of a mobile component as a 
function of pressure and temperature in a way consistent with certain assump- 
tions, In summary, these are: 

Case 1. The mobile component in question is present as a pure phase 
at the same temperature and pressure as the surrounding rock, This is the 
assumption implicit in current hydrothermal experimental work and in calcu- 
lated equilibria such as that of Goldschmidt for calcite-wollastonite. It yields 
a maximum equilibrium temperature at any given pressure (depth). 

Case 2. The rock is in chemical equilibrium with a fissure system con- 
taining a dilute aqueous solution. This method is applicable only when the 
mobile component in question is water itself. It indicates a decrease in most 
dehydration temperatures with depth. 

Case 3. The vertical gradient of chemical potential of the mobile com- 
ponent has reached a value such that there is no longer any tendency for 
vertical diffusion of this component. This method also suggests (not without 
reservations) a slight decrease of most devolatilization temperatures with 
depth. 

Case 4. The chemical potential of the mobile component is sufficiently 
low that no phases containing it are stable. 

For simplicity, we will refer to these limiting cases as cases 1, 2, 3, and 
4, respectively. It now remains for us to consider the extent to which these 
may be applicable to some specific geological problems. 

Crystallization of igneous rocks.—It is a point of general agreement 
among petrologists that the crystallization of a magma will be accompanied 
by a tendency toward supersaturation of the magma with regard to its “vola- 
tile” components. It is also a fairly safe assumption that an unaltered igneous 
rock contains less of any “volatile” component than the magma from which 
it is crystallized. This means that magma during crystallization is a “source” 
of volatiles and should tend to maintain the chemical potential of the volatile 
components at relatively high values if the “escape” of these components is 
in any way hindered. That is, we should expect relatively high temperatures 
for “devolatilization” equilibria, possibly approaching case 1 in some in- 
stances. There is, in fact, some evidence to suggest that in certain instances 
the chemical potentials of the volatiles have reached a value such that a 
“vapor” phase has separated and boiled off, and there is certainly no doubt 
that vapor phases have formed when the pressure on a magma has suddenly 
been lowered in a volcanic eruption, as many vesicular lavas bear witness. 
To the extent that a “vapor” formed in this way may be regarded as a 
one-component phase (most likely in the case of water), we have a fair ap- 
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proximation to case 1. The field of stability of biotite or hornblende, then, 
under the conditions of case 1, should not greatly exceed its field of stability 
in igneous rocks. 
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Fig. 3. (A) Possible stability relations among Al-silicates for metamorphism at 
igneous contacts (intermediate between case 1 and case 2). The field of kyanite is shown 
as masked by that of hydrous Al-silicates (pyrophyllite, kaolinite). In hydrothermal 
experiments (case 1) the fields of andalusite and sillimanite appear to be masked as 
well, Position of kyanite-andalusite-sillimanite triple point is not known, and nature of 
sillimanite-mullite transition is uncertain. 

(B) Possible stability relations among Al-silicates for regiona] metamor- 
phism (approaching case 2). 

A similar relationship among the aluminum silicate polymorphs has been 

proposed by Miyashiro (1949b). 





98 James B. Thompson, Jr—The Thermodynamic 


Metamorphism at igneous contacts—The conditions prevailing in a 
contact aureole should not differ greatly in this regard from those in the 
magma itself, inasmuch as the country rock on heating would also be a source 
of volatiles. In a limestone, furthermore, a relatively high chemical potential 
of carbon dioxide might well be obtained, though unlikely to approach case 
1, owing to the presence of magmatic water. 

Regional metamorphism.—A shale or limestone undergoing regional 
metamorphism would also constitute a source of volatiles. It would appear, 
however, that here the chances of building up or maintaining a high chemical 
potential of any given volatile would be relatively less than in the preceding 
examples. The rocks are likely to be kept “open” by orogenic movements, 
and the time involved is presumably greater. In a relatively slow metamorphic 
process, then, the mobile components have more opportunity to escape, and 
we may expect an approach to the conditions of case 2 or case 3. 

Both fissure equilibrium and free vertical diffusion indicate, in general, 
little variation of devolatization temperatures with depth, It is probably safe, 
therefore, to assume that a given devolatilization will take place at a lower 
temperature for a given depth in typical regional metamorphism than in a 
typical contact aureole, though the two may overlap inasmuch as there ap- 
pear to be all gradations between “regional” and “contact” metamorphism. 

Stress and antistress minerals——The generalization that a given “de- 
volatilization” will occur at a lower temperature, other things being equal, in 
regional metamorphism than in contact metamorphism suggests that certain 
of the “stress” minerals of Harker may be merely low temperature anhydrous, 
or nearly anhydrous phases, and that their fields of stability in “thermal” 
metamorphism may be “submerged” by those of the corresponding hydrous 
phases; that is, that the breakdown of hydrous Al-silicates to yield Al,SiO; 
may occur in the field of kyanite (fig. 3) in most regional metamorphism, 
but may be delayed until] andalusite, sillimanite or mullite is stable in contact 
metamorphism, 


A similar comment might be made with regard to the non-appearance 
of pyrope garnet in the hydrothermal experimental work of Yoder (1952) ; or 
to the observed relations between jadeite, analcite, nepheline and albite 
(Yoder, 1952; Adams, 1953) ; and possibly to the absence of an anthophyllite 
field in the work of Bowen and Tuttle (1949). In each of these instances the 
experimental work has been done under the condition Ps equals Py209 and 
should thus yield maximum dehydration temperatures. 

Metamorphism of igneous rocks—The metamorphism of pre-existing 
igneous rocks offers an interesting problem in that the initial material is 
devoid of volatiles, or nearly so, and should, if the metamorphism is at a 
sufficiently low temperature, tend to acquire them, Thus an igneous body such 
as a dunite might constitute a local “sink” for volatiles in a later metamor- 
phism. The observed peripheral alterations to serpentine and thence to talec- 
carbonate indicate that this is so, at least in some instances. Under such 
conditions, then, we may expect abnormally /ow temperatures in or near such 
bodies for equilibria involving volatiles. An igneous mass, in other words, 
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may be expected to behave as a local desiccator during subsequent low-tem- 
perature metamorphism. 

Polymetamorphism.—Polymetamorphism, at least where it involves a 
low temperature metamorphism following one at high temperature, offers a 
parallel situation to the low-temperature alteration of igneous rocks in that 
the initial material is dehydrated or partially so. Thus, in the low-grade meta- 
morphism of an older high-grade terrane, as in the regional metamorphism 
of pre-existing hornfels (Tilley, 1935), or in the Paleozoic metamorphism of 
Precambrian gneiss in the Green Mountains of Vermont, we should expect 
anhydrous mineral facies to be stable at relatively low temperatures, Noting 
that most solid-solid equilibrium curves have positive slopes, this suggests 
the possibility of the formation at abnormally low temperatures and pressures 
of facies that under the conditions of case 2 or 3 would be limited to high 
temperatures and pressures. It is tempting, in this regard, to speculate that 
some occurrences of eclogite and glaucophane schist may be of polymetamor- 
phic origin (recrystallized pyroxene hornfels) in the regions where, were 
volatiles more readily available, greenstone schists might otherwise have 
formed. Another possibility is that these are products of low temperature 
desiccation of greenstone schists inasmuch as they commonly occur in associa- 
tion with serpentinized or partly serpentinized ultramafics. 

Hydrothermal veins.—Fissure equilibrium or, in a more restricted sense, 
fissure equilibrium beneath an “impermeable barrier” may be regarded as 
an attempt to establish a model for the thermodynamic treatment of hydro- 
thermal veins and their relations to the surrounding rock, There are at least 
two restrictions, however, that might well be emphasized. One is that the 
method is most successful in dealing with reactions involving the solvent, in 
other words, to hydrations and dehydrations. To apply it to metasomatic 
processes involving other components we need more data than now available 
on the propefties of their aqueous solutions. The other restriction is that the 
temperatures of rock and adjacent fluid must not be greatly different, the 
problem otherwise being one of “steady state” conditions rather than a true 
thermodynamic equilibrium. 

Interpretation of experimental results—In a recent discussion of his 
hydrothermal investigations Yoder (1952) distinguishes between conditions 
of “excess” and “deficient” water, analogous to “oversaturation” or “under- 
saturation” with respect to silica or some other fixed component. Since Yoder 
did not consider the effect of other components in the aqueous phase, nor the 
effect of a pressure upon it less than that on the surrounding crystals, “excess 
water” would thus correspond to our case 1, and a total “deficiency” of water 
to our case 4. In a sense, case 2 and case 3 are special instances of water 
“deficiency,” but the analogy breaks down here in that we have treated rocks 
as open systems with regard to water and certain other components, whereas 
Yoder’s treatment is applicable, strictly, only to a closed system. 

Yoder, in other words, considers equilibria with regard to a given water- 
content or percent water, whereas we have done so with regard to a given 
chemical potential of water. Yoder’s variables might be given as P, T and 
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Nyeo Or % H.O and ours as P, T and py. The two, furthermore, are not 
equivalent. If a given phase containing water is present at a given temperature 
and pressure, the amount of this phase and hence the percent of water in the 
system may be increased indefinitely without altering in any way the chemical 
potential of water in the system. In another example we may have present, 
at a given pressure and temperature, a hydrate and its corresponding anhy- 
dride in equilibrium, These stipulations, alone, are sufficient to determine 
the chemical potential of water in the system, but to determine percent water 
we must know the relative amounts of hydrate and anhydride. Thus knowledge 
of P, T and number of moles of water or percent water does not define the 
chemical potential of water nor is knowledge of P, T and pxzo sufficient to 
specify the water content. 

Clearly it is an advantage to use percent water as a variable in certain 
types of experimental work, inasmuch as reactions in a “bomb” are in a 
closed system and must run at a given P (or V), T and “percent” water. In 
the Earth’s crust, however, it appears more reasonable to assume that “percent 
water” is not fixed, and that reactions will proceed at a given P, T and chemi- 
cal potential of water rather than at a given P, T and percent water. 

Experimental apparatus in which the chemical potential of water may 
be “controlled” is a possibility. One approach is to apply the hydrostatic 
pressure to the solids by means of some fluid phase other than pure or vir- 
tually pure water, thus making possible values of pyeo less than the molar 
free energy of pure water at the same pressure and temperature. This will 
lower the equilibrium temperature at a given pressure for any dehydration. If 
the “other” fluid component is an inert gas, incapable of reacting with the 
solids, its role in the experiment may be regarded as analogous to that of 
the walls of the bomb—that is, as a necessary part of the apparatus. Although 
real gases at high pressures are not “ideal” in the sense of obeying the ideal 
gas law, they are, in certain instances ideal mixtures. For water as a com- 
ponent of an ideal mixture, we have 


Luo = pH20° + RTIn Nivo (64) 


where py20° is the chemical potential of pure water at the same P and T, and 
Nuo the mole-fraction of water in the gas. If the mixture is not ideal we 
must use the relation 


Lmo = pw20° + RTin Nuvo fueo (65) 


where feo is the activity coefficient of water in the mixture. If the activity 
coefficients are known, or the system is ideal, we can thus determine px20 
from the mole-fraction of water in the fluid phase and the P-V-T data for 
pure water (Kennedy, 1950). In this regard it would be advantageous to have 
the volume of fluid large relative to that of the reacting solids so that Nu2o 
may be regarded as nearly constant. Alternatively, the equilibrium relations 
for any values of sy20 may be calculated from the experimental results with 


“excess” water if the entropies, volumes and heat capacities of the condensed 


phases are known. It is perhaps not out of place to emphasize here the definite 
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need for more thermochemical data on rock-forming minerals, and for ad- 
ditional P-V-T data on such geologically important gases as water and carbon 
dioxide. 
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RADIOACTIVITY IN GROUND AND SURFACE WATER 
SHELDON JUDSON anp J. KENNETH OSMOND 


ABSTRACT. Values for the uranium content and total radioactivity of some under- 
ground and surface waters are presented for samples from 77 localities, largely in the 
United States. The uranium content for these samples ranges between 0.02ppb and 460ppb. 
The uranium content and total radioactivity of water from carnotite-bearing beds in the 
Grand Junction, Colorado area are spectacularly high when contrasted with values for 
water from non-mineralized beds. 

Uranium content and total radioactivity was determined from residues obtained by 
evaporation of water samples. The total radioactivity of these residues is shown to change 
through time. 


INTRODUCTION 

General statement.—A portion of the geochemical cycle of uranium in- 
cludes the transportation of the element by water, both on and beneath the 
surface. Little, however, is present in the literature on the uranium content of 
such natural waters. For instance Kuroda (1953) in a compilation of data on 
the radioactivity of natural waters lists uranium values for one locality in 
Japan and two in Portugal. Furthermore, most references to the radioactivity 
of natural waters emphasize sea water, and mineral and thermal springs. For 
these waters values for radium, radon, thoron, and total activity expressed in 
radium equivalent are most common. As an exception to this, however, Adams 
(1954) has listed the uranium content of some surface waters in Wisconsin. 

The investigations here reported are primarily concerned with the uran- 
ium carried in ground water and secondarily with that in surface water. 
Brief consideration is given to the nature of other radioactive elements found 
in these waters. 

Purpose of investigation ——This study was begun to determine the amount 
of uranium carried in underground waters and, if possible, the nature of 
additional radioactive elements. The investigations were also pursued to es- 
tablish whether or not the radioactivity of ground water might be useful in 
prospecting for radioactive deposits. 

The investigation has produced some values for the uranium content of 
underground and surface waters from various environments, and some rudi- 
mentary speculations on other radio elements in these same waters. Further- 
more it seems feasible to use the radioactivity of water in geochemical pros- 
pecting under certain conditions. 

Acknowledgments.—This study is part of a general program supported 
by funds granted by the U. S. Atomic Energy Commission to the University 
of Wisconsin under University contract No. AT (11-1)-178. The project since 
its inception in 1951 has been under the general direction of Professor Far- 
rington Daniels, Chairman of the Department of Chemistry, University of 
Wisconsin. Uranium analyses in this report were made in the analytic section 
of the project under the supervision of Dr. J. A. S. Adams, now of the De- 
partment of Geology, Rice Institute. Dr. Adams has also criticized this report 
in manuscript form. The bulk of the alpha count measurements were made 
in the Department of Chemistry at Wisconsin by Mr. John Ockerman, now 
of the Ray-O-Vac Co., Madison, Wisconsin. A few of the alpha counts were 
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made by Miss Joyce Bauer, also of the Wisconsin Chemistry Department. 
Dr. M. Starr Nichols, Assistant Director of the Wisconsin State Laboratory 
of Hygiene at Madison kindly provided analyses of the fluorine content of 
certain waters. 

The writers are indebted for assistance in the collection of samples to 
Mr. William A. Oesterling, Alcoa Mining Company, and Mr. E. A. Brecke, 
Ozark-Mahoning Co., both of Roseclare, Illinois, to Dr. Henry Faul, and 
Mr. David Phoenix, both of the U. S. Geological Survey. 


METHOD OF INVESTIGATION 

Collection of samples—Samples were collected by the authors from 
various underground and surface localities in Wisconsin, Illinois, Kentucky, 
Texas, New Mexico, and the Colorado Plateau states of Utah and Colorado. 

A total of 77 water samples from various underground and surface lo- 
calities are reported on in this paper. Their general distribution is given in 
table 1. Of these 68 were 7.6-liter (2-gallon samples), 5 were 8-liter samples 
and 2 were 7-liter samples. The sample was evaporated to dryness in a steel 
pan and the residue, usually amounting to a few grams, was removed by 
scraping. It is estimated that more than 95 percent of the residue obtained 
from each sample was recovered. In all cases evaporation was carried on 
within a few days of the collection of the sample, and in many instances only 
a few hours elapsed between collection of the sample and evaporation, The 
residues were later analyzed in the laboratory. 























Taste 1 
General Distribution of Water Samples Discussed in this Report 
Location No. of Samples 

Wisconsin 50 
Southern Illinois-Western Kentucky 8 
Colorado Plateau in vicinity of Grand Junction, Colorado ............000000 8 
Texas Panhandle 6 
Canada 3 
Others 2 

Total 77 








At the time of collection the approximate pH of the water was determined 
with a colorimetric field test using phenyl red. Samples which appeared 
turbid were filtered before evaporation. 

Alpha counts.—The radioactivity of the residues was measured in terms 
of alpha counts per hour using scintillation counters.’ 

A circular sample dish 2 inches (5.08 cm) in diameter contains a 
“thick source” of the sample. The distance between the counting surface of 
the sample and the phosphor screen is 0.78 cm. The counting geometry has 
been estimated by Ockerman as about 33 percent. 

The alpha scintillation counters used were of the type described by Reed (1950). Some 
of the counters were constructed by Mr. John Ockerman at the Department of Chemistry, 


University of Wisconsin. Others were commercially manufactured units purchased from 
the Nuclear Instruments and Chemical Corp., Berkeley, California. 
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Samples were counted long enough to obtain a standard deviation of less 
than 10 percent. Most of the alpha counts reported in tables 1-4 were ob- 
tained on residues from 6 to 9 months after evaporation of the water. 

The number of alpha counts per hour of a sample is a measure of the 
radioactivity of the residue only at the time of counting. As discussed later 
the radioactivity of the sample changes with time, and it is therefore im- 
possible to obtain a check of the alpha counts of any one sample by later 
counting (table 2). Furthermore, because the sample counted was a residue 
obtained by evaporation, any radioactive gases such as thoron or radon 
which may have been present in the original water sample would not be 
present in the residue at any time, being driven off by the process of evapora- 
tion. Despite these facts the total alpha counts give at least a relative measure 
of the intensity of radioactivity in the residue and hence in the water from 
which it was derived. 

TABLE 2 
Change of Total Radioactivity of Water Residue with Time 





Sample No, 5-5-25. Eagle Hutson mine sump, Salem, Ky. 
U p.p.b. (water) = 0.5 p.p.b. Wt. residue = 2.7 g. 








Time (months) Pie yy, 9 29 
UG, BUD. crecccnesecnsencnconesonsescoseesnsnesutinenesvanatesnestnes 34 13 15 
CEPR ODD cctsnsevsecreasshithenstniuasceostisibtetaninngiielictlilndiis 1% 


Sample No, 9-3-1, Cactus Rat Mine seep, Colorado Plateau 
U p.p.b. (water) = 100 p.p.b. Wt. residue = 11.4 ¢. 
Time (months) 2 6 23 
U-eq. p.p.b. 370 520 380 
U/U-eq. (%) ... 27 19 26 











Uranium content of the residue.—The uranium content of the individual 
samples of residue was determined under the supervision of Dr. Adams using 
the fluormetric method described by Grimaldi, May and Fletcher (1952). 
Repeated analyses were found to fall within the range of a standard error 
originally estimated to be between 10 and 15 percent. 

Uranium content of the water—The uranium content of the original 
water sample was obtained by the following formula: 


Uppb Uppm X gms. residue x 10° 
(water ) (residue) gms. water evaporated 


Total radioactivity of the water.—Uranium accounts for only a portion 
of the radioactivity of the original water. An attempt was made to determine 
the total radioactivity of the water in terms of uranium equivalent using the 
following formula: 


U-equiv. ppb alpha cph X2Xx gms. residue C 103 
{water ) (residue) ~ "gms. water evaporated — 


The factor of “2” was determined empirically in the laboratory where 
it was found that standard water samples gave one alpha cph for 2 ppm 
of uranium. 

The values of uranium-equivalent, however, are at best an approxima- 
tion of the total radioactivity of the original water. In addition to the errors 
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inherent in residue recovery in alpha counting and in uranium recovery other 
considerations 'imit the precision of the final figures of uranium-equivalent. 
Thus radioactivity of the residue changes through time as already mentioned, 
and thus the alpha counts represent only a measure of the radioactivity of 
the residue at time of counting. Secondly, radioactive gases in the original 
water would be lost in evaporation and cannot be represented in the resulting 
residue. 

Therefore the values of uranium-equivalent given for the original water 
sample are approximations only. The same is true of the percentage of the 
radioactivity of the original water due to uranium. This value is obtained by 
dividing the uranium content of the water by the uranium equivalent of the 
radioactivity of the water in ppb and multiplying by 100. Nevertheless the 
figures probably have a qualitative value although their degree of precision 
is difficult to judge. 

l'rocedural checks and conclusions.—Procedural checks were made in 
the laboratory with standard and “spiked” water samples to determine (1) 
the proportion of residue recovered after evaporation of the water, (2) the 
proportion of the uranium recovered from the residue, and (3) the alpha 
counts produced by a given quantity of uranium in the residue, These checks 
showed (1) that at least 95 percent of the residue is recovered by scraping, 
(2) that at least 70 percent of the uranium is recovered from the residue 
sample, and (3) that 2 ppm of uranium yields one alpha count per hour. 

On the basis of these facts and in view of the fair consistency of the 
values for the uranium content for different samples of waters of the same 
general type, it is felt that the uranium content of residue and original water 
as given in the tables to follow have quantitative significance. The validity 
of these values is further emphasized by the fact that Adams (1954) obtained 
values for the uranium content of Wisconsin surface waters similar to those 
reported in table 9. Adams obtained his data by the fluormetric analysis of 
water samples rather than by the analysis of residues as in this study. 

The values in uranium equivalent for total radioactivity of the original 
water and the values for the percentage of uranium of total radioactivity in 
the original water are appproximations only, as previously indicated. Despite 
this these values are probably of the correct order of magnitude. 

The data in the following sections are reported with the full appreciation 
of the limitation imposed by the method of the investigation here reported. 


URANIUM CONTENT AND RADIOACTIVITY OF WATER 


Data for 73 of the 77 water samples collected are listed in tables 3 through 
9. The four samples not tabulated are discussed in the text. 

Ground water from the Colorado Plateau.—Table 3 lists data on waters 
from aquifers in the uranium mining district of the Colorado Plateau in the 
vicinity of Grand Junction, Colorado. Five samples were collected from 
mineralized (carnotite) aquifers and three from non-mineralized aquifers. 
The table clearly shows that the waters from the mineralized areas are more 
radioactive than those from the non-mineralized aquifers. In uranium con- 
tent waters from mineralized aquifers range from 100 to 460 ppb, whereas 








TABLE 3 


Data on Samples from Localities 


Sheldon Judson and J. Kenneth Osmond 


on the Colorado Plateau in Colorado and Utah 
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| | = 
= a 
| | an ee 
— ————_—_—_—_—_—__——— OO Se i a 2 
Sam | & 3 z Se] &ix 
ple | 4 ae px a Mt ta 
No. Location Source Aquifer pH ei & S168 S| slit 
ate | | @l|2lesi | eie 
‘s s |e /8°] sj sic 
Mineralized Aquifer s 5 z xsi >| DID 
9.3-1- Cactus Rat Mine, mine “Morrison 8.3 114 66 173 6 100 520 19 
Thompson, Utah seep fm. 
9-4-1- Yellow Cat Camp, spring Morrison >84 4.1 468 664 6 250 720 35 
Thompson, Utah fm. 
9-5-1 Yellow Circle Camp, spring Morrison >8.4 25 640 454 6 200 300 67 
Moab, Utah fm. 
9-8-1 45-90 Mine, mine Morrison >84 3.7 950 627 6 460 610 75 
Uravan, Colorado seep fm. 
9.9.1 Henry Clay mine, mine Morrison >84 2.0 1040 1767 6 270 930 29 
Uravan, Colorado sump fm. 
Non-mineralized Aquifer 
9-6-1 Joe Davis Canyon, spring Morrison ~ >84 74° 21 , §} =~ at 2 
Egnar, Colorado fm. 
9-7-1 Gypsum Gap, spring Dakota >84 94 7.2 SB 6 98 @ LU 
Naturita, Colorado fm. 
9-10-1 Dolores River spring Wingate >84 19 109 So 6 27. © 2 


Canyon, Gateway, 


fm. 


Colorado 


in non-mineralized aquifers the range is from 2 to 9 ppb. The total radioac- 
tivity and the percentage of radioactivity due to uranium are also significantly 
higher in the mineralized than in the non-mineralized zones. All waters were 
alkaline having pH values of 8.3 or more. 

Ground water from northeastern Wisconsin.—Ground-water samples from 
Cambrian sandstone in northeastern Wisconsin are listed in table 4, and 
water from Ordovician, Silurian and late Pleistocene aquifers of the same 
area are listed in table 5. 

The waters from the Cambrian sandstone were obtained from wells which 
were drilled far into the aquifer and which, in some cases, reached to the 
Precambrian basement. The Paleozoic rocks of the area dip to the southeast 
off the Wisconsin arch or dome of Precambrian rocks. Therefore the wells 
increase in depth toward the southeast and bottom between 200 feet and 1000 
feet beneath the surface. 

Compared with the Colorado Plateau waters, particularly those from 
mineralized aquifers, these Wisconsin waters are much lower in uranium 
content, ranging from 0.1 to 1.7 ppb. The total radioactivity expressed as 
uranium-equivalent is much lower than that of the water from the carnotite- 
bearing beds of the Colorado Plateau but is similar to the total radioactivity 
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Data on Samples from Deep Wells in Northeastern Wisconsin 


Aquifer in all instances, Cambrian sandstone. 
Well depths range between 200 and 1000 feet and average about 800 feet. 





Residue 














Uranium ppm 

Alpha cph 

Time elapsed 

to count (months) 
Uranium ppb 
Uranium equiv. ppb 
U/U-eq. (%) 


Sample 
0. 


Location Well 




















10-1-1 Appleton Outagamie Asylum 
10-2-1 Kaukauna City Well #4 
10-3-1 Wrightstown Fox River Dairy 
10-8-1 Point Sauble Village Well 
10-9-1 Pulaski City Well 
10-10-1 Suamico Suamico Dairy 
10-13-1 Bonduel City Well 
10-14-1 Seymour City Well #2 
10-17-1 DePere East Bank Well 
10-18-1 De Pere Armour & Co. 
10-19-1 Green Bay Ashwaubenon 

Sanitary District 
10-20-1 Green Bay Gray Street Well : 
10-21-1 Green Bay Boland Road Well . 0 1.91 
10-22-1 Preble Deckner Ave. Well ; 8 2.69 . 
10-23-1 Preble Town Well #2 i 4 5.54 7 130 
10-24-1 Green Bay Northern Paper Mills 3 2.89 a 3 
10-25-1 Hortonville City Well F 1 2.38 ’ 36 
10-26-1 Marinette Marinette Paper Mills 2 5.44 é 93 
10-27-1 Peshtigo Peshtigo Paper Mills .7 1.43 9 119 
10-29-1 Lena Village Well } 1 1.43 7 102 
10-30-1 Oconto Falls City Well #3 y .1 1.44 8 91 
10-31-1 Oconto City Well #6 . 6 1.67 109 
10-32-1 Coleman Village Well , .1 2.00 i 67 
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of water from non-mineralized beds of that area. The total radioactivity of 
the water due to uranium is low when compared to Colorado Plateau waters 
from both mineralized and non-mineralized aquifers. 

The pH of the waters from the Wisconsin Cambrian rocks ranges between 
7.4 and 8.2. There is some indication that the total amount of uranium 
carried in these waters decreases with increasing pH. In this regard we may 
note that Phair and Levine (1953) report that uranium in the form of UO, 
is preferentially leached in H,SO, solutions at the expense of UO., Ra and 
Pb and further that increasing concentration of H.SO, increases the propor- 
tion of UO, leached. 

The ground water in northeastern Wisconsin is relatively high in fluorine 
ranging from 0.1 ppm to 2.8 ppm. The correlation between uranium and 
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TABLE 5 
Data on Samples from Wells of Intermediate Depth in 
Northeastern Wisconsin. Aquifers as Indicated 
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Sample . Eis & -| < + 
No. Location Well pH | = a ee ae a ~ a 
Aquifer: Ordovician ss. Depth: 200 to 1300 feet; Average, 500 feet 
10-41 Denmark Chi. & N.W. RR Depot 84 814 43 4 6 46 9 50 
10-5-1 _ Brillion City Well 76 215 18 10 6 05 6 8 
10-6-1 Algoma City Well #2 714 382 2441 6 12 4! 3 
10-11-1 Suamico Carl Jenkins 8.2 166 90 42 6 20 18 11 
10-12-1 Mill Center Walter Kimbs 8.4 087 03 6 6 03 14 2 
10-15-1 Kaukauna Rapid Crocke Powerhouse 7.6 6.30 08 54 6 07 9% 1 
10-16-1 Little Rapids Brown Co. Sanitarium 8.0 280 05 18 6 0.2— iM a reee 
Aquifer: Ordovician limestone. Depth: 30 to 150 feet. “" 
10-37-1 De Pere Harry Liebergen ~ 84 166 04 6 6 O1 2% 4 
10-38-1 Little Rapids Roy Wellmer 82 3.19 08 6 6 03 5 6 
10-47-1 Seymour George Sieserman 74 236 06 4 6 O02 2% 8 
Aquifer: Silurian Is. Depth: 70 to 106 feet. : 
10-7-1 Casco Junction Green Bay and 7.4 256 46 35 6 16 24 6 
Western RR Well 
10-35-1 Algoma Louis Manly 18 2% 15 6 6 UG 4 13 
10-36-1 Dykesville John Collin 718 198 13 + © GA 3% 11 
10-39-1 Shirley George Treml 76 7.98 15 16 6 16 33 5 
Aquifer: Glacial Drift. Depth: 50 to 100 feet. 
10-40-1 Green Bay Clarence Greatens 82 3.06 01 3 6 005 2% 2 
10-42-1 Pound Leo Cudnoski 84 159 65 24 6 14 10 14 
10-43-1 Oconto Falls Walter Birr 7.6 209 06 8 6 O02 4% 4 
10-45-1 Embarrass Anthony Anton 8.4 0.75 04 11 6 005 2 2% 
Aquifer: fractured granite and glacial drift. Depth: 350 feet. , 
10-46-1 Embarrass Wedig Krubsach 8.0 142 100 29 6 19 Ill 17 














fluorine is not marked, but in general the total radioactivity of the water in 
uranium-equivalent increases with an increase in the fluorine content of the 
water, The fluorine concentration is greatest in the deeper wells down the 
dip of the Cambrian sandstone with the highest values recorded along the 
Fox River Valley between Lake Winnebago and Green Bay. 


The results of water measurements from shallower and stratigraphically 


higher aquifers than the Cambrian in northeastern Wisconsin are given in 
table 5. The aquifers are Ordovician sandstone and dolomite, Silurian dolo- 
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mite, and Pleistocene (Wisconsin) gravel. The waters from these aquifers are 
similar in uranium content to that from the Cambrian, but their total radio- 
activity is not so great. The proportion of radioactivity due to uranium is 
variable but averages about 5 or 10 percent. 

Ground water from the Texas High Plains —Table 6 lists the data on 
five water samples from the Ogallalah formation (Pliocene) of the High 
Plains of Texas. Their uranium content, 2 to 10 ppb, is intermediate between 
that of the Wisconsin waters and the waters from mineralized waters of the 
Colorado Plateau. They are similar in uranium content to the waters from 
the non-mineralized aquifers of the Colorado Plateau. Their total radioac- 
tivity is small but the percentage of activity due to uranium is between 25 
and 63, similar to the water from the carnotite-bearing aquifers of Colorado 
and Utah. The source of uranium in these samples is thought to be volcanic 
ash which occurs in Pleistocene beds above the Ogallalah and probably in 
the Ogallalah itself. Ground-water recharge of the Ogallalah is local. 


TABLE 6 
Data on Samples from the Texas Panhandle 


Aquifer: Ogallalah fm. (Pliocene) 


Residue 
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No. Location Source 


Uranium ppm 
Alpha eph 
Uranium ppb 
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8-3-3 Exell, Texas Stock well — il 19 
8-4-1 Exell, Texas Stock well ; ‘ 16 
8-7-1 Dawn, Texas House well : 25 44 
8-7-3 Hereford, Texas Irrigation well . 21 20 
8-8-3 Tule Canyon, Texas House well 03 25 20 
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Ground water from southern Illinois-western Kentucky.—The results for 
seven water samples from the southern Illinois and western Kentucky Fluospar 
district are listed in table 7. Six of these samples were collected from mine 
workings, and one from an artesian spring. The uranium content is uniformly 
low, having values between 0.13 and 0.5 ppb. 

Miscellaneous ground water samples.—A sample from the warm spring 
at Pagosa Springs, Colorado, and a sample from an artesian well at Newkirk, 
New Mexico, were collected. The results are given in table 8. 

Surface waters——Waters from surface streams were collected at nine 
localities, seven of them in Wisconsin. Data on these waters are given in table 
9. The uranium content is generally low as is the total radioactivity in terms 
of uranium equivalent. The percentage of radioactivity due to uranium tends 
to be high. 

Waters not tabulated.—The data for four additional samples are not 
included in tabular form. One of these was collected from a spring at the 
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TABLE 7 
Data on Samples from Southern Illinois and Western Kentucky 
Aquifer: Mississippian Limestone 





Residue Water 





Sample 
No. Location Source 


Uranium ppm 
Alpha cph 

Time elapsed 

to count (months) 
Uranium ppb 


























| Uranium equiv. ppb 
>| U/U-eq. (%) 


5-3-2 Rosiclare, Ill. Rosiclare Fluor- 
spar Mine 
Cave-in-Rock, Ill. W. Green Fluor- 
spar Mine 
Cave-in-Rock, Ill. Deardorf Fluor- 
spar Mine 
Rosiclare, Ill. Alcoa Fluor- 
spar Mine 
Salem, Kentucky Eagle Hutson 
Zine Mine 
Cave-in-Rock, Ill. “Rd. Spring” 
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* 8 liter sample. 


Fission Mines, Wilberforce, Ontario; another from a pond in a pegmatite 
quarry, Madawaska, Ontario. The third, an 8-liter sample, was supplied by 
Mr. Alan Gregory from Beaverlodge Lake in the Lake Athabaska Area, 
Saskatchewan. Each of the samples was chosen because of proximity to a 
uraniferous source. In each instance, however, the amount of residue re- 
covered was so small that the analytical results were not judged to be reliable. 
A fourth sample of water from Cambrian sandstone at Madison, Wisconsin 
is referred to later. 


TABLE 8 
Data on Miscellaneous Samples 





Residue 





Sample 
No. Location Source ph 


8-9-1 Newkirk,N.M. Housewell >8.4 9.88 


9-1-1 Pagosa Springs, (Triassic ss) 7.0 19.10 
Alo. Warm Spring 
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Conclusions.—On the basis of the data presented above the following 
conclusions can be drawn: 

1. Both underground water and surface water carry uranium in measur- 
able amounts. It cannot be said whether all of this uranium is in solution or 
whether it is in suspension associated in some way with colloidal-size particles. 

2. Where underground waters pass through an aquifer containing carn- 
otite, the total radioactivity and uranium content are high when compared 
with water from similar but non-mineralized aquifers, It cannot be said on the 
basis of the data collected whether the same relationships would hold for 
aquifers containing uranium-bearing minerals other than carnotite. On the 
other hand, the relatively high uranium content and radioactivity of waters 
from carnotite-contaminated beds suggests that geochemical methods might 
well be devised for prospecting for uraniferous deposits. 

3. The residue technique used in this study appears to be practical only 
for those waters which carry much material in solution and thus produce a 
relatively large residue sample, say 1 gram per 2-gallon water sample. 

4. The relation between pH and total uranium carried is not clear from 
the data at hand. Waters from Cambrian rocks in northeastern Wisconsin, 
however, decrease in uranium content as the pH increases from 7.4 to 8.0. 

5. Incomplete data obtained from water in Cambrian rocks of north- 
eastern Wisconsin also suggest that in these waters fluorine increases with 
increasing radioactivity but not necessarily with uranium. 


CHANGES IN RADIOACTIVITY WITH TIME 
It has already been stated that the alpha counts per hour for any one 
residue sample change with time. The effect of this change on the calculations 
of total radioactivity of the water in uranium-equivalent and of the percentage 
of radioactivity due to uranium has been illustrated in table 2. Successive 
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Fig. 1. Decrease in radioactivity with time thought to be due to thorium X (radium 
224 in thorium series). Sample is a residue of water from the Cambrian sandstone at 
Madison, Wisconsin, taken at the main station well. The sample counted at successive 
intervals, Large scatter of points due largely to short counting periods. 
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recounts of a single sample showed that the change of radioactivity was 
greatest and least predictable in the first few months after evaporation. 
Samples counted 20 or more months after evaporation were found to be 
changing at a slower rate. 

An excellent example of decreasing radioactivity in terms of alpha counts 
per hour is provided by a sample of residue obtained from the evaporation 
of water from a deep well in Cambrian sandstone at Madison, Wisconsin. 
The sample was counted repeatedly during the three-week period immediately 
following evaporation of the water. The activity decreased by one-half during 
this interval as shown in figure 1. 

Table 10 lists 25 water samples for which the residue of each was counted 
two or more times at varying intervals. The data are not complete enough to 
warrant any definite conclusions or generalizations except (1) that the 
changes are real, (2) the radioactivity in terms of alpha counts may either 
increase or decrease with time, and (3) the direction of change appears to 
reverse itself in some instances. 

Members of the 4n + 1 series cannot be cited as responsible for the 
change because they are not known to occur in nature, The actinium series 
is exceedingly rare in nature so must be excluded from consideration. This 
leaves members of the uranium and thorium series. Of these, elements with 
long half lives such as Uranium 238 and 234 and Thorium 232 and 230 
disintegrate at too slow a rate to be reflected over the short time span of a 
few weeks, or months, or at most something over two years during which 


TABLE 9 
Data on Surface Waters 
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No. Location Source pH 


Alpha cph 

Time elapsed 

to count (mont 
Uranium ppb 
Uranium equiv. ppb 























5-7-3 Cave-in-Rock, Il]. Deardorf Creek 7.2 3.09 
5-2-4 Rosiclare, Ill. Ohio River 7.8 2.65 
8-6-1 N. of Amarillo, Canadian River >8.4 4.67 
Texas 
10-28-1 Middle Inlet, Wis. Granite Quarry >8.4 1.50 
Seep* 
10-33-1 Mountain, Wis. Oconto River >84 262 38 
10-34-1 Clintonville, Wis. Pigeon River >8.4 1.41 16.0 
10-44-1 Shawano, Wis. Wolf River >8.4 2.64 1.2 
10-41-1 Green Bay, Wis. Fox River 78 122 29 
10-48-1 Appleton, Wis. Fox River >8.4 1.06 3.7 
10-49-1 Madison, Wis. Lake Mendota >8.4 2.86 1.1 
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Taste 10 
Change in Radioactivity of Water Residue with Time in Terms 


of Alpha Counts per Hour 





Time Time Time Time 
elapsed elapsed elapsed elapsed 
tocount Alpha| tocount Alpha] tocount Alpha] to count Alpha 

No. Location (months) cph |(months) _cph (months) cph |(months) cph 
5-3-2 Rosiclare, Il. 7 10 9 10 
5-4-1 Cave-in-Rock, Ill. 7 121 9 68 18 48 
5-5-25 Salem, Ky. Y%y 51 9 19 29 22 
5-6-9 Cave-in-Rock, Il. y, 5 9 4 
9 
9 





5-7-1 Cave-in-Rock, III. \, 20 10 
5-8-1 Rosiclare, Il. 13 
8-3-3 Exell, Texas 0 98 
8-4-1 Exell, Texas 0 56 
8-7-1 Dawn, Texas 0 50 
8-7-3 Hereford, Texas 
8-8-3 Tule Canyon, Texas 
8-9-1 Newkirk, N. M. 

9-1-1 Pagosa Springs, Colo. 
9.3-1 Thompson, Utah 
9-4-1 Thompson, Utah 
9.5-1 Moab, Utah 

9-8-1 Uravan, Colo. 

9-9-1 Uravan, Colo. 

9-6-1 Egnar, Colo. 

9-7-1 Naturita, Colo. 
9-10-1 Gateway, Colo. 
10-10-1 Suamico, Wis. 
10-14-1 Seymour, Wis. 
10-18-1 De Pere, Wis. 
10-21-1 Green Bay, Wis. 
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627 
1767 
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successive counts were made. The cause of the change then must lie among 
the members of the uranium and thorium series with shorter half lives. These 
would include those elements with half lives of Thorium X (radium 224) or 
less, but exclusive of radon and thoron which were not present in the original 
residue sample. Furthermore, we may exclude those elements of half lives 
less than three or four days because of the impracticability of measuring them 
by the alpha count method. 

A series of activity curves for members of the uranium and thorium 
series was constructed. Theoretically the comparison of these curves with the 
changes in alpha activity of the water residues through time should permit 
speculations concerning the relative importance of some of the daughter ele- 
ments of uranium and thorium. Unfortunately, data are insufficient for most 
samples to allow any precise plotting of the nature of the change of alpha 
activity through time. Only in the case of a sample obtained from waters in 
the Cambrian sandstone at Madison, Wisconsin are enough alpha determina- 
tions available to sketch the change in alpha activity through time. This 
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sample was counted at successive intervals over a three-week period and the 
results given in figure 1. It is presumed that the rapid decrease in activity 
during the first two weeks was due primarily to the presence of thorium X 
(radium 224 in the thorium series). 
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NOTES ON THE ORIGIN OF 
SOME AGATES AND THEIR BEARING ON 
A STYLOLITE SEAM IN PETRIFIED WOOD 


B. M. SHAUB 


ABSTRACT. The “channels,” common in agate, that have been considered to be influx 
canals are believed to be the result of outward moving bulbous masses of the more gelat- 
inous and plastic silica gel after rhythmic banding has developed in the segregated 
material. A similar plastic gelatinous state appears to exist during an intermediate period 
of development of at least some petrified wood. During this gelatinous stage stylolites or 
stylolite-like structures can develop in the still plastic material under certain conditions 
where masses of the plastic gel can move into zones of reduced pressure formed during 
and as a result of the dehydration and consolidation that yielded the jasper or chalcedony 
of the petrified wood. 


The theories of origin of agates have been discussed on numerous oc- 
casions, and it is beyond the scope of this paper to review these here. There 
appear, however, to be a few features that have been overlooked or misin- 
terpreted, and as these have a bearing on the nature of a piece of petrified 
wood in which a stylolite seam has developed, it is appropriate to discuss 
briefly the probable origin of agates like the one illustrated in figure 1. 

In 1915 Liesegang described many agates as natural examples of a 
process of rhythmic diffusion after the manner in which he had produced 
rhythmic banding in gelatinous preparations by the diffusion of reagents 
causing colored precipitates. In these preparations the reacting agents entered 
the gelatinous preparation from the exterior. Liesegang presumed that in the 
banding of agates the coloring and banding material entered the agates in a 
similar manner and circulated through certain tubes, canals and other pass- 
ages in the silica gel. 

Configurations such as those shown at a, b, c, d, and e, figure 1, were 
described by J. Noeggerath (cited by Liesegang, p. 17) as “einflusskaniile” 
which were supposed to be influx canals through which material entered the 
agate. To the writer, such an operation appears to be contrary to the facts, 
as shown by the configuration of the agate bands associated with these “ca- 
nals.” These show that the so-called influx canals are in reality exit canals 
instead, The bulbous material at a, b, c, d, and e has moved outward to its 
present position from zones near the interior of the mass. 

The sharply banded and fine-grained chalcedony is succeeded in the 
interior of the agate nodule by a crust of poorly terminated quartz crystals. 
The structure and texture of the agate suggest the following conclusions: (1) 
the banding occurred in a mass of gelatinous material; (2) the banding 
developed prior to the movement of the material from near the interior of 
the mass to a position near the periphery;: (3) the crust of quartz crystals 
was probably the last to form although some of the gelatinous and nearly 
transparent silica remained in the center. 

Agates of this type occur chiefly in the fine-grained felsitic or porphyritic 
gabbroic or more basic rocks as nodules varying in size up to 10 or 12 inches 
but usually smaller. Judging from the characteristics of the original material, 
as shown by the features of the agate illustrated, the gelatinous silica mass 
must have separated as an immiscible fraction in the basic magma toward 
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Fig. 1. Agate showing exit channels and contiguous bulbous masses of chalcedony 
at a, b, c, d, e and elsewhere, which have moved outward from the interior of the agate 
structure while the entire mass was still soft. 


the end of its consolidation. It probably remained as a gelatinous mass in- 
closed in the solidified rock for some time before the banding, and other 
features were developed, as deduced from this and numerous other specimens. 

In such material it is not difficult to find a logical motive force to cause 
the transfer of material from the interior to the peripheral parts of the mass, 
for in such a gelatinous mass inclosed in a now granular or partly granular 
rock the water in the exterior portions of the gelatinous mass would slowly 
escape or be dispersed in the surrounding rock, and the outside material 
would contract an amount equal to the volume of water removed, The more 
aqueous interior portions of the mass, now under slightly higher pressure, 
would then slowly move outward along the zones or channels of least resist- 
ance, producing the bulbous parts and simultaneously bending and dragging 
the previously formed bands outward in the direction of the moving material 
and, in so doing, reducing the thickness of the bands by the dragging action 
in the direction of motion. 
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In this agate at the present time we have a material that, except for its 
banding, is extremely close in its composition to the present material in petri- 
fied wood both as to mineralogical and physical properties. 

The above discussion has been presented to indicate that certain agates 
originally consisted of a plastic gelatinous material capable of being deformed 
by flowage. 

In the case of petrified wood too, we have numerous theories of origin 
which cannot be discussed here. It has appeared to the author in discussing 
the petrification of wood with other geologists that the assumed processes 
which yield the silicified wood, according to any of the theories, would pro- 
duce at first a solid chalcedonic mass that is hard and rigid. Accounts of 
recently silicified wood or of wood in the process of silicification appear to 
be rare; this also reflects the scarcity of such material available for observa- 
union. 

Reporting of field work in Patagonia, Argentina, Nichols and Miller 
(1951, p. 282) observed that a “piece of wood that had been altered to a 
transparent jelly was found in the weathered eolian silt between the forest 
litter and the layer of volcanic ash. It is not known whether this was the 
result of the normal decay of buried wood under the conditions existing in 
this part of the world or whether the wood had been replaced by silica gel 
derived from the ash and is therefore an intermediate stage in the formation 
of petrified wood.” 


By 


Fig. 2. A polished surface on a section of petrified wood showing an irregularly 
developed stylolite seam along the periphery or outer portion. The bulges at a, b, and 
c are essentially expanded “columns” while the long narrow one at d extends deeply into 
the silicified wood. The monochrome photograph loses many of the details of the pene- 
trating parts which are more clearly shown by their color differentiations. 











120 B. M, Shaub 


It is to be noted that in this instance the partly silicified wood was pass- 
ing through the stage of a translucent jelly capable of being deformed under 
low stresses either from contraction within the gelatinous mass by gradual 
loss of water to the surrounding materials or by direct external forces, 

A few years ago a piece of polished petrified wood was brought to my 
attention that contained around the outer edge an irregular stylolite or stylo- 
lite-like seam in which there are a number of unusual bulgy “columns” (a, 
b, and c, fig. 2), also a long projecting “column,” d, of chalcedony penetrating 
deeply into the specimen. That these projections extending into the mass of 
the petrified wood could have originated according to the pressure-solution 
theory for the origin of stylolites, or by any of the more recently modified 
forms of this theory, is indeed unacceptable; first, on account of the highly 
insoluble nature of the hardened material; second, there was never enough 
thickness of solid chalcedony about the bark periphery to provide columns 
of the length shown; and third, the bulgy character of the columns in their 
hardened condition would not permit movement of the irregular, jagged and 
interlocked parts in any direction. Although the polished section is cut ap- 
proximately normal to the axis of the tree trunk or limb, it is not likely that 
the movement of the material of the bulges or the long column, d, was exactly 
centripetal in nature but that it moved by plastic flow in many directions 
depending upon the release of the pressure in the surrounding material. Each 
section, across or at various angles to the axis of the tree, would be expected 
to show a unique pattern either with or without bulged columns as shown in 
the particular section illustrated. 

Such stylolitic structures, however, can originate according to the con- 
traction-pressure theory in silicified woods during the intermediate gelatinous 
phases of their development when the loss of water would cause a shrinkage 
and contemporaneous plastic flow of adjacent more yielding silica gel into 
the zones of reduced pressure. 
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SPECTROGRAPHIC DETERMINATION OF 
“ONTAMINATION OF ROCK SAMPLES 
APica, ot [NDING WITH ALUMINA CERAMIC* 


?, R. BARNETT, W. P. HULEATT, L. F. RADER, ann A. T. MYERS 


ABSTRACT. Contamination of rock samples pulverized between experimental ceramic 
disks is determined and found to be almost negligible when compared to contamination 
resulting when the samples are ground with steel disks, particularly with regard to iron 
and alloying elements of steel. A method of mounting the disks on a conventional grinder 
is described. 

A previous study by Myers and Barnett (1953) showed that rock samples 
pulverized in heavy grinding machinery were contaminated with tramp iron 
containing certain alloying elements used in the manufacture of steel. A 
search was started for suitable materials to substitute for steel grinding plates. 
Arrangements were made with the Coors Porcelain Company, Golden, Colo- 
rado, to make experimental grinding plates from a high-alumina base; this 
material is relatively pure, particularly with regard to alloying elements of 
steel, and is hard and durable. It is probable that high-alumina ceramic 
products of other manufacturers would be equally satisfactory. One set of 
alumina-ceramic grinding plates ground about 1400 two-ounce rock samples, 
mostly sandstones, to -80 mesh before replacement was necessary. Thus the 
grinding of rock samples with a ceramic material was successfully demon- 
strated as a means of rapid preparation of samples for spectrographic analysis 
and is now standard practice in one grinding laboratory of the U. S. Geologi- 
cal Survey. 

Figure 1 shows the grinder with a set of plates attached and ready for 
use. The bottom plate is attached to the vertical rotating shaft of the grinder 
by means of a tapered steel bushing through the center of the plate. This plate 
is supported and strengthened by an aluminum housing machined to hold it 
firmly in place. Two steel plugs with threaded holes are inserted diagonally 
opposite each other in the upper plate near its edge to receive bolts that hold 
the plate stationary when mounted in the grinder. A second set of ceramic 
grinding plates is shown in the right foreground of figure 1. The pan to catch 
the ground sample is also shown. This pan is made of aluminum to prevent 
sample contamination by other metals. 

The kind and degree of contamination introduced by the ceramic plates 
was determined by grinding selected massive vein quartz and quartzite samples 
from the same source as the samples used in the previous study (Myers and 
Barnett, 1953). As in the previous study, control portions of the rocks were 
ground to —100 mesh in agate and analyzed spectrographically to compare 
with splits pulverized to -100 mesh with the ceramic grinding plates. The 
method of analyses was that described by Fleischer and others (1952). 

The quantitative data in table 1 show the effectiveness of the ceramic in 
preventing contamination of rock samples, especially with respect to iron and 
the other elements usually found in steel. 

*Publication authorized by the Director, U. S. Geological Survey. 


121 





‘saqejd JIWIBIIO YIM Surpuns 49]jB JUBUIMIBJUOD BSB Pojoejop JOU PUB [O1}UOD ay) Ul peeep you JN 10j peyoo, 949M yeud s}yuvgul 
“3/9 4940 40} os|e ‘sqIWwTy] AVA ISUaS 40} zZ 37qe) 99S *pevelep you nq “104 pe yoo, S}UIUII9 — 0 +10} peyoo] you SyUIUII]o == *apqey aaoqe uy 
‘auryoeu Surpurss yi jo yyeys 9atup oy) UO Surysng 19ddoo ay) 07 paoesy SBM UOTJBUTWIEUOD Jaddoo Stuy, » » 


*(€S61) NeuIeg pue sidkyy Wo weg , 


. 
™— 
~ 
5 
= 
~ 
~ 
S 
— 
° 
~ 
s 
~ 
— 
~ 
= 
= 
_ 
~ 
~ 
= 
~ 
v 
~ 
. 
ms) 
— 
2 
= 
<= 
~ 
— 
~ 
~ 
_ 
SO 
= 
L 
~ 
— 
* 
—— 
“,) 





L000°0 
900°0 
100 
c00°0 
a | 
100°0 
100°0 
£000°0 


# [9995 


0 
0 
0 
0 


« « LOOO'O 


0 
0 


£000°0 
0 

0 
£00°0 
0 


0 
LOOO'O 
0 
0 


JTUIBIAGD YIM 


ayede ul 
: Surpus 


£000°0 
£00°0 
900°0 
L00°O 
60 
£000°0 
6000°0 
Z000°0 


. jeer 





0 
0 
0 
0 
0 


0 
0 
0 
0 
0 
0 
0 
« «£000°0 
0 
0 


oTUIB189 


10°0 
£00°0 
600 
600°0 
b0°0 
10°0 
c0°0 
F0°0 


F000°0 
0 
9000°0 
0 

Ol 
£000°0 
&000°0 
0 


OTUIBI9D YIM 


100 
S00°0 
b0'0 
c00°0 
600 
200°0 
600 

600 
£0000 

0 

9000°0 

0 

OT 

0 

£000°0 

0 

oayese ul 
:Suipuns 


=z 
A 
ML 


Jaye uotjeuruajap ‘oadg saye uoneurmajep ‘ved¢ 


:Aq UOTyBUTUIE UO) 


:Aq uoljeuTWEUO) , 


zuenb daisseyy ayzyien?) 


(juaosad ul) jaaIg YIM pue o1uresa7) YIM SuIpuL4 Aq Zeng) OAISSE]Y 


P.R. Barnett. et al. 


pue ayizyen¢) JO uoTeUTUIR]UO’) Surseduio’) sisdjeuy orydesZo1jo0dg aarjeyuen?) 
[ FMavy 


122 





nation of Rock Samples after Grinding with Alumina Ceramic 123 


— 0.02 
— 0,001 
- 0.001 
0.0001 
— 0.0001 
- 0.0005 
0.0005 
0.001 
0.005 
— 0.0001 
— 0.0001 
- 0.00005 
0.05 
0.005 
0.001 
0.0001 
- 0.005 
0.0005 


TABLE 2 
Sensitivities of Elements Analyzed for by the 
Techniques of this Study 


Hf 
Hg 
In 
Ir 
La 
Li 
Mg 
Mn 


Mo - 


Na 
Nb 


Nd — 


Ni 
Os 
P 
Pb 
Pd 
Pt 
Re 


% 


f —0.05 — 


— 1.0 
— 0.001 
— 0.005 
— 0.003 
— 0.01 
— 0.001 
~— 0.0001 
- 0.0001 
— 0.02 
0.0005 
0.01 

— 0.0001 
— 0.005 





Grinder with set of ceramic plates attached. A second set of plates is in the foreground. 
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This work is part of a program conducted by the U. S. Geological Survey 
on behalf of the Division of Raw Materials of the U. S. Atomic Energy 
Commission. 

The authors are designing a new type of grinder for use with the alumina- 
ceramic grinding plates for pulverizing rock samples under conditions that 


are even less subject to metal contamination than the present system. 
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REVIEWS 


Principles of Geomorphology; by WittiaM D. THornsury. P. ix, 618; 
274 figs. New York, 1954 (John Wiley & Sons, $8.00).—This timely book 
fills a genuine need in American college teaching today, which is for a text 
in geomorphology written for the upper division geology major just starting 
his professional courses..-The book does an excellent job of meeting precisely 
that need, It evidently presupposes the education and experience of a junior 
or senior student who has had one or more beginning courses in geology. 

The specialist unquestionably will find points to quarrel over with the 
author, and undoubtedly will feel pain over the brevity of treatment afforded 
his particular field of interest—if so, the book is not .vr him. Professor Thorn- 
bury has done outstandingly well in steering a middle course—he has not 
taken the direction of writing at too elementary a level, nor has he veered in 
the direction of attempting to impress a small, select group of his peers with 
his erudition. 

Geomorphology now includes so many specialties that it is a courageous 
man indeed who attempts to encompass the entire field. On the whole, the 
author does succeed in this endeavor. Necessarily the treatment of individual 
topics is thin, but this is inevitable in a work of this sort. From my point of 
view it would have been desirable had more extensive treatment been given 
some of the controversial topics of the day, such as submarine canyons, coral 
reefs, the origin of pediments, and the nature of slope retreat. 

To compensate for the brevity of the discussion of many branches of the 
science, the book includes a wide-ranging and unusually comprehensive list 
of references at the end of each chapter. These should prove a godsend to 
the interested student of the subject, as well as to the disinterested one if he 
is required to write a term report. The references have been carefully selected, 
and this is especially true for controversial issues, in order to bring out more 
than one point of view. The references are up to date, and were they to be 
read, would give a person a good understanding of the present status of 
geomorphology. The specialist undoubtedly will deplore the omission of a 
favorite article, especially if he is the author, but this is to be expected in a 
work such as this. The standards used for inclusion seem to have been almost 
completely unbiased. 

In fact, it is the lack of bias that makes this such a satisfactory text. 
This is not to say that students should be or are sheltered from conflict, but 
in the case of a recognized controversy, the various positions are stated fairly, 
and then instead of leaving the student with an unsolved mystery on his hands, 
the explanation most generally accepted today is clearly set forth. 

The historic growth of the science, and the contribution of individuals 
to that growth are well handled. I wish personally that even more emphasis 
had been placed on this approach because geomorphology is a branch of geol- 
ogy whose present-day aspect is the work of a relatively small number of 
remarkably strong-minded and highly individualistic men. It should be in- 
teresting to the student to see how the science has developed as a result of the 
interplay of their ideas. 

The photographic illustrations are an excellent addition to the text. For 
the most part they are aerial obliques; the majority are new, and all go far 
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in showing the true appearance and mutual relationships of landforms, The 
line drawings, while adequate, are not up to the standard we have grown to 
expect from such gifted men as Cotton, Davis, King, Lobeck, and Raisz— 
but this is the sort of achievement few of us are in a position to duplicate. 

An interesting chapter on Applied Geomorphology closes the book. In 
it are cited examples taken from economic geology, engineering geology, 
petroleum geology, and military geology in which a knowledge of geomorph- 
ology has proved useful. The examples used by no means encompass the entire 
field, and it would be valuable indeed to see this chapter expanded in later 
editions. 

The book is well made and is priced within the limits of the average 
student’s budget. It suffers from its share of typographical errors, few of 
which are not immediately obvious, and one involving a chapter heading will 
be mildly amusing to all but the author. 

This book is an excellent introduction to the field of geomorphology, 
and with the wealth of references that are included—most of which are avail- 
able in the average departmental library—there is every opportunity to lead 
the student into an enduring interest in the subject. 


WILLIAM C, PUTNAM 


The Permian Reef Complex of the Guadelupe Mountains Region, Texas 
and New Mexico: A Study in Paleoecology; by Norman D. NEWELL, J. KEITH 


Ricsy, ALrrep G. Fiscner, A. J. WairemMan, Joun E, Hickox, and JOHN 
S. Brapiey. P. xix, 236; 85 figs., 32 pls. San Francisco, 1953 (W. H. Freeman 
& Company, $7.50).—The Permian reefs of West Texas have been intensively 
studied for decades, and we thought we knew a great deal about them. The 
recent study by Newell and his colleagues has added much more, until now 
the growth of the Capitan reef is one of the best known chapters in geologic 
history older than the Ice Age, yet they think of their book only as a progress 
report and are the first to point out how much still remains to be learned. 
Their findings on these ancient reefs are the more valuable because several 
members of the team prepared themselves by studying modern reefs in the 
Bahamas and the Pacific. 

The theme of their book is best indicated by its subtitle: A Study in 
Paleoecology. There are short supporting chapters on the regional setting, 
stratigraphy, and paleontology, but nearly half the book is a long chapter 
describing in detail the eleven phases or facies of the reef complex, from the 
sandstone in the center of the basin through the reef proper to the sandstone 
on the landward edge of the lagoon. Though these descriptions are more 
petrologic than ecologic, the intimate interrelation of rock facies and life 
facies is stressed throughout, both the control of organisms by their environ- 
ments, and the creation of new environments by the reef-building organisms. 
This reviewer was particularly interested in the description of the reef-talus 
phase, which is difficult of access in modern reefs and hence has received 
relatively little attention in most recent studies of modern atolls and barriers. 


Yet it forms half or more of the massive Capitan formation, the core of the 
reef complex. 
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Another important chapter is that on diagenesis. Here the emphasis is 
on dolomitization and especially silicification—both silicified fossils and chert 
nodules. Convincing explanations of the origin of each are built up from 
fundamental chemical principles. Final chapters summarize the paleoecology 
of the reef complex and discuss where within it oil originates and accumulates. 

In an early footnote (p. 11), the authors reject polysyllabic terminology 
on principle, and indeed the whole book demonstrates very well how clearly 
complex scientific material can be presented in plain English, if only sufficient 
care be taken, (That writing in polysyllables can be hopelessly obscure re- 
quires, unhappily, no demonstration.) There are however a few blemishes in 
the writing; for example, the authors have occasionally followed a popular 
trend among geologists and have hidden the verb- idea of a sentence in a 
noun ending i in -tion (“Distinction is readily made. . .”). 

No review of the book would be complete that did not heartily thank 
the Humble Oil and Refining Company for their farsighted and objective 
support of this important project, and W. H. Freeman & Company for publish- 
ing in such excellent form (note the 32 collotype plates) a book that would 
seem at first sight to have a very limited market. In fact, the book is so full 
of valuable information about reefs and limestones generally that any geologist 
who deals with sedimentary or at least carbonate rocks must have a copy. The 
two companies have made possible a major contribution to the science of 
geology; may they also reap their financial reward. 

JOHN RODGERS 


Neuere Probleme der Abstammungslehre. Die transspezifische Evolution. 
Zweite, stark veranderte Auflage; by BERNHARD Renscu. P. xi, 436; 113 
figs. Stuttgart, 1954 (Ferdinand Enke Verlag, DM 47.-, paper cover; 49.20, 
clothbound) .—The first edition of this book appeared in 1947. At that time 
the author did not have access to the extensive foreign and especially Ameri- 
can literature on the same subject, and for this reason his book did not receive 
much attention outside of Germany. The second edition is greatly superior 
in this respect. Several chapters have been completely rewritten, especially 
those dealing with the tempo and mode of evolution and the one concerning 
the phenomena of consciousness. Many of the 113 figures illustrating the text 
have been borrowed from American publications and are familiar to Ameri- 
can readers, A few of the illustrations are original. Twenty-six pages of 
closely printed references to published books and articles bear evidence of 
the great erudition of the author and of his thorough acquaintance with the 
main subject of “transspecific” evolution. This term may be somewhat mis- 
leading. Rensch objects to the term macroevolution as being an awkward 
combination of Greek and Latin and has coined the term “transspecific” for 
the evolution of all categories which are higher than the species. To the re- 
viewer Rensch’s term seems to be even more unfortunate than the older term, 
inasmuch as it could be equally well if not better applied to subspecific evolu- 
tion, whereas, “macroevolution” is free of such confusion, In referring to 
“the evolution of the small and large branches of the main stem” Rensch uses 
the term “cladogenesis” and proposes a new term “anagenesis” for “higher 
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development” as exemplified by the ever-growing complexity of structure in 
all branches of the phylogenetic tree, He makes extensive use of special terms 
and seems sometimes to forget that their value depends entirely upon the 
extent to which they can help us to understand the phenomena for which they 
are coined, and that their frequent use, while giving a book the appearance 
of profound erudition, makes the text less readable and thus limits the circle 
of possible readers. 

As a zoologist the author stresses primarily the evidence derived from 
the study of those branches of biological sciences which deal with living 
animals. His approach to the problems of evolution will appeal to those who 
believe that speciation is the only mode and the true basis of all evolution. 
Those who like the present reviewer think that macroevolution cannot be 
explained by the process of speciation, will find the book deficient in the dis- 
cussion of material which should have been included. Thus Rensch gives 
barely three pages to the question of sexual isolation and, although he himself 
states that “purely ecological isolation may lead to speciation,” the section on 
ecological races gives very little exact information. Being an adherent of the 
idea that all mutations responsible for the production of the great multiplicity 
of animal forms are in themselves non-directional, the author explains the 
direction of evolution as a result of natural selection. Yet in the case of ana- 
genesis he is not satisfied with this explanation and expresses the opinion 
that higher development is of necessity directional. 

The last chapter of the book is devoted to a discussion of the evolution 
of psychic phenomena. The author raises four questions which he tries to 
answer. These questions are: (1) Is it probable that an evolution of the phe- 
nomena of conciousness has taken place? (2) How far down the animal scale 
can the existence of sensations, perceptions, mental processes and voluntary 
actions be traced? (3) Can somatic evolution be influenced by psychic phe- 
nomena? and (4) How does the evolution of the latter fit into a general 
philosophical concept of the world? It seems evident that no definite answer 
could be given to these questions. Perhaps the author would have been wiser 
not to raise them at all, as an adequate survey of animal behavior and of the 
problems of psychic phenomena would in itself require more than a single 
volume. In every other respect the present book is a mine of information. It 
is accompanied by a good subject index and should be of real use to a student 
of the problems of evolution. 

ALEXANDER PETRUNKEVITCH 








Longwell and Flint'’s new geology . . . 


INTRODUCTION to 
PHYSICAL GEOLOGY 


By Chester R. Longwell, Henry Barnard Davis Professor of Geology, 
and Richard Foster Flint, Professor of Geology, both of Yale University. 


A streamlined, authoritative presentation . . . 

This modern introduction concentrates on the principles 
of physical geology and develops them through a thorough 
examination of actual field evidence. Evidence bearing on 
the age of the Earth’s crust is presented early and often, 
as are the concepts of igneous geology. The formation of 
fresh “primary” rocks is logically considered in advance 
of the changes wrought by weathering and erosion, The 
conflict of forces responsible for the “perpetual motion” 
of erosion, sedimentation, and uplift of lands is emphasized 
repeatedly in the book, which aims at setting the geologic 
processes in true perspective for the beginning student. 

Clear and graphic... . 

Representatives of more than 140 college geology de- 
partments contributed ideas to the arrangement and relative 
emphasis of subjects in this book. It is distinguished by its 
use of simple language and only a minimum of technical 
terms, More than 340 significant photographs and drawings 
accompany the text, many with printed labels and explana- 
tory diagrams, Critical parts of all drawings are labeled 
directly, and graphic scales are provided wherever dimen- 
sions are important. These features illuminate an introduc- 
tory geology that is concise, clear, and economical. Send for 
an on-approval copy. 
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